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° SIGNAL PROCESSING APPARATUS 

Backcrround of the Invention 
Reference to Prior Related App lication 

This is a continuation-in-part application of U.S. 
5 Patent Application No. 08/132,812 filed October 6, 1993, and 

entitled "Signal Processing Apparatus" and a continuation-in- 
part application of U.S. Patent Application No. 08/249,69-0 
filed May 26, 1994 entitled "Signal Procesfi-ing-A$>pa-r-a-tu«-and 
Method" which is a continuation of U.S. Patent Application 
10 No. 07/666,060 filed March 7, 1991, now abandoned. 
Field of the Invention 

The present invention relates to the field of signal - 
J processing. More specifically, the present invention relst&a 

to the processing of measured signals, containing a primary 
I 15 signal portion and a secondary signal portion, for the 

I removal or derivation of either the piluidiy oir secuadaiy 

i signal portion when little is known about either of these 

components. More particularly, the present invention relates 
I to modeling the measured signals in a novel way which 

20 facilitates minimizing the correlation between the primary 

: signal portion and the secondary signal portion in order to 

1 produce a primary and/or secondary signal . The. nrpspnt 

invention is especially useful for physiological monitoring- 
systems including blood oxygen saturation systems. 
25 Description of t he Related Art 

Signal processors are typically employed to. remove or 
derive either the primairy or secondary signal portion from a 
composite measured signal including a primary signal portion 
and a secondary signal portion. For example, a composite 
30 signal may contain noise and desirable portions. If the 

secondary signal portion occupies a different frequency 
spectrum than the primary signal portion, then conventional 
filtering techniques such as low pass, band pass, and high 
pass filtering are available to remove or derive either the 
35 primary or the secondary signal portion from the total 
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signal. Fixed single or multiple notch filters could also be 
employed if the primary and/or secondary signal portion (s) 
exist at a fixed frequency (s) . 

It is often the case that an overlap in frequency 
5 spectrum between the primary and secondary signal portions 
exists* Complicating matters further, the statistical 
properties of one or both of the primary and secondary ^signal- 
portions change with time. In such cases^, conventional 
filtering techniques are ineffective in extracting eit liex Lire " 

10 primary or secondary signal. If, however, a description of 
either the primary or secondary signal portion can 
derived, correlation canceling, such as adaptive rrnrrere-" 
canceling, can be employed to remove either the primary or 
secondary signal portion of the signal isolating the other 

15 portion. In other words^ given sufficient information -about - 

one of the signal portions, that signal portion can be 
extracted. 

Conventional correlation cancelers, such as adaptive 
noise cancelers, dynamically change their transfer function 

20 to adapt to and remove portions of a composite signal. 

However, correlation cancelers require either a secondary 
reference or a primary reference which correlates to either 
the secondary signal portion only or the primary signal 
portion only. For instance, for a measured signal containing 

25 noise and desirable signal, the noise can be removed with a 

correlation canceler if a noise reference is available. This 
is often the case. Although the amplitude of the reference 
signals are not necessarily the same as the amplitude of the 
corresponding primary or secondary signal portions, they have 

30 a frequency spectrum which is similar to that of the primary 

or secondary signal portions. 

In many cases, nothing or very little is known about the 
secondary and/or primary signal portions. One area where 
measured signals comprising a primary signal portion and a 

3 5 secondary signal portion about which no information can 

easily be determined is physiological monitoring. 
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Physiological monitoring generally involves measured signals 
derived from a physiological system, such as the human body. 
Measurements which are typically taken with physiological 
monitoring systems include electrocardiographs, blood 
5 pressure, blood gas saturation {such as oxygen saturation), 
capnographs, other blood constituent monitoring, heart rate, 
respiration rate, electro- encephalograph (EEG) and depth of 
anesthesia, for example. Other types of measurements include 
those which measure the pressure and quantity of a substance 

10 within the body such as cardiac output, venous oxygen^ 
saturation, arterial oxygen saturation, bilirubin, total 
hemoglobin, breathalyzer testing, drug testing, cholesterol 
testing, glucose testing, extra vasation, and carbon dioxide 
testing, protein testing^ carbon monoxide testing, anfi ofhfir 

15 in-vivo measurements, for example. Complications arising in 

these measurements are often due to motion of the patient, 
both external and internal (muscle movement, vessel movement, 
and probe movement, for example), during the measurement 
process. 

20 Many types of physiological measurements can be made by 

using the known properties of energy attenuation as a 
selected form of energy passes through a medium, 

A blood gas monitor is one example of a physiological 
monitoring system which is based upon the measurement of 

25 energy attenuated by biological tissues or substances. Blood 

gas monitors transmit light into the test medium and measure 
the attenuation of the light as a function of time. The 
output signal of a blood gas monitor which is sensitive to. 
the arterial blood flow contains a component which is a 

30 waveform representative of the patient's arterial pulse. 

This type of signal, which contains a component related to 
the patient's pulse, is called a plethysmographic wave, and 
is shown in Figure 1 as curve s. Plethysmographic waveforms 
are used in blood gas saturation measurements. As the heart 

3 5 beats, the amount of blood in the arteries increases and 
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decreases, causing increases and decreases in energy 
attenuation, illustrated by the cyclic wave s in Figure 1. 

Typically, a digit such as a finger, an ear lobe, or 
other portion of the body where blood flows close to the 
5 skin, is employed as the medium through which light energy 

is transmitted for blood gas attenuation measurements ♦ The 
finger " comprises skin, fat, bone, muscle, etc., shown 
schematically in Figure 2, each of which att:eircs:t:s^ eii^di'^y - 
incident on the finger in a generally predictable and 

10 constant manner* However, when fleshy portions of the f inaf»- 

are compressed erratically, for example by motion of the 
finger, energy attenuation becomes erratic. 

An example of a more realistic measured waveform S i"s 
shown in Figure 3, illustrating the effect of motion. The 

15 primary plethysmographic waveform pQr^^on of. thf^ signal s.-ls 

the waveform representative of the pulse, corresponding xro- 
the sawtooth-like pattern wave in Figure 1, The large, 
secondary mot ion -induced excursions in signal amplitude 
obscure the primary plethysmographic signal s. Even small 

20 variations in amplitude make it difficult to distinguish the 

primary signal component s in the presence of a secondary 
signal component n, 

A pulse oximeter is a type of blood gas monitor which 
non-invasively measures the arterial saturation of oxygen in 

25 the blood. The pumping of the heart forces freshly 

oxygenated blood into the arteries causing greater energy 
attenuation. As well understood in the art, the arterxal 
saturation of oxygenated blood may be determined from the 
depth of the valleys relative to the peaks of two 

3 0 plethysmographic waveforms measured at separatre wavel-e-ii y ths '. 

Patient movement introduces motion artifacts to the composit-er 
signal as illustrated in the plethysmographic waveform 
illustrated in Figure 3. These motion artifacts distort the 
measured signal. 
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Summary of the Invention 
This invention provides improvements upon the methods 
and apparatus disclosed in U.S. Patent Application 
No, 08/132,812, filed October 6, 1993, entitled Signal 
5 Processing Apparatus, which earlier application has been 

assigned to the assignee of the instant application. The 
present invention involves several different embodiments 
using the novel signal model in accordance wi t h the pzsjb ' e::L • 
invention to isolate either a primary signal portion or a 

10 secondary signal portion of a composite measured signal. Ixl 
one embodiment, a signal processor acquires a first measured 
signal and a second measxired signal that is correlated to the 
first measured signal. The first signal comprises a first 
primary signal portion and a first secondary signal portion. 

IS The second signal comprises a second primary signal port i on 

and a second secondary signal portion. "The signals may be 
acquired by propagating energy through a medium and measuring 
an attenuated signal after transmission or reflection. 
Alternatively, the signals may be acquired by measuring 

20 energy generated by the medium. 

In one embodiment, the first and second measured signals 
are processed to generate a secondary reference -which does 
not contain the primary signal portions from either of the 
first or second measured signals. This secondary reference 

25 is correlated to the secondary signal portion of each of the 

first and second measured signals. The secondary reference 
is used to remove the secondary portion of each of the first- 
and second measured signals via a correlation canceler, such 
as an adaptive noise canceler. The correlation canceler is 

3 0 a device which takes a first and second input nnd removes 

from the first input all signal components which are- 
correlated to the second input. Any unit which performs or 
nearly performs this function is herein considered to be a 
correlation canceler. 

35 An adaptive correlation canceler can be described by 

analogy to a dynamic multiple notch filter which dynamically 
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changes its transfer function in response to a reference 
signal and the measured signals to remove frequencies from 
the measured signals that are also present in the reference 
signal. Thus, a typical adaptive correlation canceler 
5 receives the signal from which it is desired to remove a 

component and receives a reference signal of the undesired 
portion. The output of the correlation canceler ±s a. good, 
approximatioa to the desired signal with the i::iid:efi±r=ad- 
component removed. 

10 Alternatively/ the first and second measured signals mav 

be processed to generate a primary reference which does not 
contain the secondary signal portions from either of the 
first or second measured signals* The primary reference may 
then be used to remove the primary portion o'f~ each of" the 

15 first and second measured signals via a correlation canceler • 

The o utput "of — the" coiielaLlun canuelex ±s^ *a gooU " 

approximation to the secondary signal with the primary signal 
removed and may be used for subsequent processing in the same 
instrument or an auxiliary instrument. In this capacity, the 

20 approximation to the secondary signal may be used as- a 

reference signal for input to a second correlation caiiceler- 
together with either the first or second measured signals for 
computation of, respectively, either the first or second 
primary signal portions. 

25 Physiological monitors can benefit from signal 

processors of the present invention. Often in physiological 
measurements a first signal comprising a first primary 
portion and a first secondary portion and a ^second signal 
comprising a second primary portion and a second secondary 

30 portion are acquired. The signals may be acquired ixy 

propagating energy through a patient's body (or a material 
which is derived from the body, such as breath, blood, or 
tissue, for example) or inside a vessel and measuring an 
attenuated signal after transmission or reflection^ 

35 Alternatively, the signal may be acquired by measuring energy 

generated by a patient's body, such as in 
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electrocardiography. The signals are processed via the 
signal processor of the present invention to acquire either 
a secondary reference or a primary reference which is input 
to a correlation canceler, such as an adaptive noise 
5 canceler* 

One physiological monitoring apparatus which benefits 
from the present invention is a monitoring system which 
determines a signal which is representative of the arterial 
pulse, called a plethysmographic wave. This, sig nal can be. 

10 used in blood pressure calculations, blood constitueiiL 

measurements, etc, A specific example of such a use is in 
pulse oximetry. Pulse oximetry involves determining the 
saturation of oxygen in the blood. In this configuration, 
the primary p ortion of the signal is the- arte-rial bleed 

15 contribution to attenuation of energy as it passes through a 

portion of the hndv-^whare ..h3oo,d..£l.ows close to the skin. The 
pumping of the heart causes blood flow to increase and 
decrease in the arteries in a periodic fashion, causing 
periodic attenuation wherein the periodic waveform is the 

20 plethysmographic waveform representative of the arterial 

pulse. The secondary portion is noise. In accordance with 
the present invention, the measured signals are modeled such 
that this secondary portion of the signal is related to the 
venous blood contribution to attenuation of energy as it 

25 passes through the body. The secondary portion also includes 

artifacts due to patient movement which causes the venous 
blood to flow in an unpredictable manner, causing, 
unpredictable attenuation and , corrupting * the oth^rwi-c^- 
per iodic plethysmographic waveform. Respi rat i on al so caiiaes- 

30 the secondary or noise portion to vary, although typically -at 

a lower frequency than the patients pulse rate. Accordingly/ 
the measured signal which forms' a plethysnaographic waveform, 
is modeled in accordance with the present invention such that 
the primary portion of the signal is representative of 

3 5 arterial blood contribution to attenuation and the secondary 

portion is due to several other parameters. 



-7- 



A physiological monitor particularly adapted to pulse 
oximetry oxygen saturation measurement comprises two light 
emitting diodes (LED's) which emit light at different 
wavelengths to produce first and second signals, A detector 
5 registers the attenuation of the two different energy signals 

after each passes through an absorptive media, for example a 
digit such as a finger, or an earlobe* The attenuated 
signals generally comprise both primary (arterial attenuator) 
and secondary (noise) signal portions. A static filtering 

10 system, such as a bandpass filter, removes a portion of the- 

secondary signal which is outside of a known bandwidth of" 
interest, leaving an erratic or random secondary signal 
portion, often caused by motion and often difficult to* 
remove, along with, the primary signal port i on - 

15 A processor in accordance with one embodiment of che 

present invention removes the primary signal portions from 
the measured signals yielding a secondary reference which is 
a combination of the remaining secondary signal portions. 
The secondary reference is correlated to both of the 

20 secondary signal portions- The secondary reference and at 

least one of the measured signals are input to a correlation 
canceler, such as an adaptive noise canceler, which xeiuoves — 
the random or erratic portion of the secondary signal. This 
yields a good approximation to a primary plethysmographic 

25 signal as measured at one of the measured signal wavelengths. 

As is known in the art, quantitative measurements of the 
amount of oxygenated arterial blood in the body can be 
determined from the plethysmographic signal in a variety of 
ways , 

30 The processor of the present invention may also remove 

the secondary signal portions from the measured signals 
yielding a primary reference which is a combination of the 
rentaining primary signal portions. The primary reference is 
correlated to both of the primary signal portions. The 

35 primary reference and at least one of the measured signals 

are input to a correlation canceler which removes the primary 
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portions of the measured signals. This yields a good 
approximation to the secondary signal at one of the measured 
signal wavelengths. This signal may be useful for removing 
secondary signals from an auxiliary instrument as well as 
5 determining venous blood oxygen saturation. 

In accordance with the signal model of the present 
invention, the two measured signals each having primary and 
secondary signal portions can be related by coefficients. By 
relating , the two equations with respect to coefficients 

10 defined in accordance with the present invention, the 
coefficients provide information about the arterial oxygen 
saturation and about the noise (the venous oxygen saturation 
and other parameters) . In accordance with this aspect of the 
present invention, the cx:>ef f icients can be df^^erm^npd by 

15 minimizing the correlation between the primary and secondary- 

signal portions as defined in the model. Accordingly, the 
signal model of the present invention can be utilized in many 
ways in order to obtain information about the measured 
signals as will be further apparent in the detailed 

20 description of the preferred embodiments. 

One aspect of the present invention is a method for use 
in a signal processor in a signal processor for piucesslua a tr 
least two measured signals and each containing a primary 
signal portion s and a secondary signal portion n, the 

25 signals and being in accordance with the. following 

relationship: 

5, = n, 

where s^ and Sj, and n^ and nj are related by: 

and where r^ and are coefficients. 

The method comprises a number of steps, A value of 
3 0 coefficient r, is determined which minimize correlation 

between s^ and n^ . Then, at least one of the first and second 

-9- 



signals is processed using the determined value for to 
significantly reduce n from at least one of the first or 
second measured signal to form a clean signal. 

In one embodiment, the clean signal is displayed on a 
5 display. In another embodiment/ wherein the first and second 

signals are physiological signals, the method further 
comprises the step of processing the clean signal to 
determine a physioTogical parameter from the first or second 
measured signals. In one embodiment, the jjarameter is 

10 arterial oxygen saturation. In another embodiment, t h^ 

parameter is an ECG signal. In yet anoLhez " euibjjdlmenL, ' 
wherein the firsu portion of the measured signals ±5- 
indicative of a hear t plethysmograph, the method further^ 
comprises the, f^tep of caTnulating the, nnlgg^ -rahA. 

15 Another aspect of the present invention invoiv-es a- 

physiological monitor. The monitor has a first input 
configured to receive a first measured signal having a 
primary portion, s^, and a secondairy portion n^^ The monitor 
also has a second input configured to received a second 

20 measured signal S2 having a primary portion S2 and a secondary 

portion nj. Advantageously, the first and the second 
measured signals and are irt accordance — wi-th~ tire 
following relationship: 

- + n^ 

where and s^, and n^ and n2 are related by: 

Si = r^S2 and n^^r^^ ' ' 

25 and where r^^ and r^ are coefficients. 

The monitor further has a scan reference processor, the 
scan reference processor responds to a plurality-xrf-'pcrssl-ble 
values for r^ to multiply the second measured signal by each 
of the possible values for r^ and for each of the resulting 

30 values, to subtract the resulting values from the first 

measured signal to provide a plurality of output signals. A 



corre lation canceler having a first input configured to 
receive the first measured signal, and having a second input 
configured to receive the plurality of output signals from 
the saturation scan reference processor, provides a plurality 
5 of output vectors corresponding to the correlation 
cancellation between the plurality of output signals and the 
first measured signal. .An integrator having an input 
configured to receive the plurality of output vectors from 
the correlation canceler is responsive to the plurality cf 

10 output vectors to determine a corresponding power tor- esch- 
output vector. An extremum detector is coupled at its input 
to the output of the integrator. The extremum detector is 
responsive to the corresponding power for each output vector 
to detect a s elected power. 

15 In one embodiment, the plurality of possible values 

correspond tCL_a_^iliirality q£ possible valnfts t or a sp.l,ect£d_ 
blood constituent. In one embodiment the, the selected blood 
constituent is arterial blood oxygen saturation. In another 
embodiment, the selected blood constituent is venous blood 

20 oxygen saturation. In yet another embodiment, the selected 

blood constituent is carbon monoxide - 

Another aspect of the present invention involves a 
physiological monitor. The monitor has a first input 
configured to receive a first measured signal having a 

25 primary portion, s^, and a secondary portion, n^. The monitor 

also has a second input configured to received a second 
measured signal Sj having a primary portion Sj and a secondary 
portion nj . The first and the second measured si-grrals -Si-artd 
S2 are in accordance with the following relationship: 

= + 

30 where s^ and Sj, and n^ and n^ are related by: 

Sj^ = r^s^ and n^^rji^ 



and where r, and r^ are coefficients. 
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A transform module is responsive to the first and the 
second measured signals and responsive to a plurality of 
possible values for r^ to provide at least one power curve as 
an output. An extremum calculation module is responsive to 
5 the at least one power curve to select a value for r^ which 

minimizes the correlation between s and n, and to calculate 
from the value for r^ a corresponding saturation value as an 
output, A display modxile is responsive to the output .of. 
saturation calculation to display the saturation valu-e-. 
10 "Brief Description of the Drawincrs 

Figure 1 illustrates an ideal plethysmographic waveform. 

Figure 2 schematically illustrates a typical finger. 

Figure 3 illustrates a plethysmographic waveform which 
includes a mouiuii - iuduced erratic signal pnrtiDrr, 
15 Figure 4a illustrates a schematic diagram of a 

physiological monitor to con^^ute priraari' — phy-si^olxiigicaJL. 
signals. 

Figure 4b illustrates a schematic diagram of a 
physiological monitor to compute secondary signals. 
20 Figure 5a illustrates an example of an adaptive noise 

canceler which could be employed in a physiological monitor, 
to compute primary physiological signals. 

Figure 5b illustrates an example of an adaptive noise 
canceler which could be employed in a physiological monitor, 
25 to compute secondary motion artifact signals. 

Figure 5c illustrates the transfer function of -a- 
multiple notch filter. 

Figure 6a illustrates a schematic of abso/:biag lUcLLeilctl 
comprising N constituents within the absorbing material, 
30 Figure 6b illustrates another schematic of absorbing 

material comprising N constituents, including one mixed 
layer, within the absorbing material. 

Figure 6c illustrates another schematic of absorbing 
material comprising N constituents, including two mixed 
35 layers, within the absorbing material. 
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Figure 7a illustrates a schematic diagram of a monitor, 
to compute primary and secondary signals in accordance with 
one aspect of the present invention. 

Figure 7b illustrates the ideal correlation canceler 
energy or power output as a function of the signal 
coefficients r^, r^, ••• r„. In this particular example, r^ 
= r. and r7 = 

Figure 7c illustrates the non- ideal correlation cancelen 
energy or power output as a function of the- signal 
coefficients r^, r,, r„. In this particular example, -ry 

= r, and r, = rv 

Figure 8 is a schematic model of a joint process 
estimator comprising a least-squares lattice predictor and a 
regression filter. 
15 Figure 8a is a schematic model of a joint process 

estimator comprising- -a QRD least-squares lattice -{LSI.). 
predictor and a regression filter. 

Figure 9 is a flowchart representing a subroutine for 
implementing in software a joint process estimator as modeled 

20 in Figure 8. 

Figure 9a is a flowchart representing a subroutine for 
implementing in software a joint process estimator as modeled 
in Figure 8a. 

Figure 10 is a schematic model of a joint process 
25 estimator with a least-squares lattice predictor and two 
regression filters. 

Figure 10a is a schematic model of a joint pro cess 
estimator with a QRD least-squares lattice pre'dictor Hxrd-tvrcr 

regression filters. 

Figure 11 is an example of a physiological monitor in 
accordance with the teachings of one aspect of the present 
invention. 

Figure 11a illustrates an example of a low noise emitter 
current driver with accompanying digital to analog converter. 
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Figure 12 illustrates the front end analog signal 
conditioning circuitry and the analog to digital conversion 
circuitry of the physiological monitor of Figure 11. 

Figure 13 illustrates further detail of the digital 
5 signal processing circuitry of Figure 11. 

Figure 14 illustrates additional detail of the 
operations performed by the digital signal proceyying" 
circuitry of Figure 11. 

Figure IS illustrates additional detail regarding the 
10 demodulation module of Figure 14. 

Figure 16 illustrates additional detail regarding tire- 
decimation module of Figure 14. 

Figure 17 represents a more detailed block diagram of^ 
the operations of the statistics module of Figure 14 . 
15 Figure 18 illustrates a block diagram of the operations 

of one embodiment of the saturation transform module of 
Figure 14 . 

Figure 19 illustrates a block diagram of the operation 
of the saturation calculation module of Figure 14 . 
20 Figure 20 illustrates a block diagram of the operations 

of the pulse rate calculation module of Figure 14 . 

Figure 21 illustrates a block diagram of the operations - 
of the motion artifact suppression module of Figure 20, 

Figure 21a illustrates an alternative block diagram for 
25 the operations of the motion artifact suppression module of 

Figure 20. 

Figure 22 illustrates a saturation transform curve in 
accordance with the principles of the present, inventioiu 

Figure 23 illustrates a block diagram of an alternative 
30 embodiment to the saturation transform in order to obtain a 

saturation value. 

Figure 24 illustrates a histogram saturation transform 
in accordance with the alternative embodiment of. Figure 23. 
Figures 25A-25C illustrate yet another alternative 
35 embodiment in order to obtain the saturation. 
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Figure 26 illustrates a signal measured at a red 
wavelength Xa=Xred=660nni for use in a processor of the 
present invention for determining the secondary reference 
n' (t) or the primary reference s' (t) and for use in a 
5 correlation canceler. The measured signal comprises a primary 

portion Sj^(t) and a secondary portion nx.(t) . 

Figure 27 illustrates a signal measured at an infrared 
wavelength Xb=Xiii=910nm for use in a processor of the present 
invention for determining the secondary reference n' (t) ox 
10 the primary reference s' (t) and for use in a correlation 
canceler. The measured signal comprises a primary portion 
SxbCt) and a secondary portion nxbCt) . 

Figure 28 illustrates the secondary reference n' (t) 
determined by a processor of the present inventiors- 
15 Figure 29 illustrates a good approximation s"xa{t) to the- 

primary portion Sj.^{t) of the signal Sj^(t) measured at 
Xa=Xred=660nm estimated by correlation cancellation with a 
secondary reference n' (t) . 

Figure 30 illustrates a good approximation s"xb{t) to 
20 the primary portion s^{t) of the signal Sxb(t) measured .at 

Xb=XIR=910nm estimated by correlation cancellation with a 
secondary reference n' (t) . 

Figure 31 depicts a set of 3 concentric electrodes, 
i.e., a tripolar electrode sensor, to derive 
25 electrocardiography (ECG) signals, denoted as Sj, Sj and Sj, 
for use with the present invention. Each of the ECG signals 
contains a primary portion and a secondary portion. 

Detailed Description of the Invention 
The present invention involves a system which utilizes 
30 first and second measured signals that each contain a primary 

signal portion and a secondary signal portion. In other 
words, given a first and second composite signals 
Si(t) = Si{t) + ni(t) and S^Ct) = s^Ct) + n^Ct), the system of 
the present invention can be used to • isolate either the 
35 primary signal portion s(t) or the secondary signal portion 

n(t). Following processing, the output of the system 
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provides a good approximation n" (t) to the secondary signal 
portion n(t) or a good approximation s"(t) to the primary 
signal portion s{t). 

The system of the present invention is particularly 
5 useful where the primary and/or secondary signal portion n(t) 

may contain one or more of a constant portion, a predictable 
portion, an erratic portion, a random portion, etc. The 
primary signal -approximation (t) or -secondary signal 
approximation n" (t) is derived by removing as many of the 

10 secondary signal portions n(t) or primary signal portions 
s{t) from the composite signal S{t) as possibi-e-. The- 
remaining signal forms either the primary sxgnal 
approximation s" (t) or secondary signal approxiraalrian n" 
respectively. The constant portion and predictable portion 

15 of the secondary signal n{t) are .easily removed v±th- 

traditional filtering techniques, such as simple subcraction, 
low pass, band pass, and high pass filtering. The erratic 
portion is more difficult to remove due to its unpredictable 
nature. If something is known about the erratic signal, even 

20 statistically, it could be removed, at least partially, from 
the measured signal via traditional filtering techniques. 
However, often no information is known about the erratic 
portion of the secondary signal n(t). In this case, 
traditional filtering techniques are usually insufficient. 

25 In order to remove the secondary signal n(t), a signal 

model in accordance with the present invention is defined as 
follows for the first and second measured signals Si and Sj: 
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5, = s,+xi, 



with 



= r^Sj and n^ = n^ 



_ ^1 



= -i; and = — 
* s. ^ n. 



where and are at least somewhat (preferably 

substantially) uncorrelated and Sj and are at least 
somewhat (preferably substantially) uncorrelated. The tirst 
and second nmcxbu.Lcia - signals and S2 r&lart-ed - by 

correlation coefficients r, and r, as defined above. The use 
and selection -of these coeffirienrs is ripsrrib ed in further, 
detail below. 

In accordance with one aspect of the present invention, 
this signal model is used in combination with a correlation 
canceler, such as an adaptive noise canceler, to remove .or 
derive the erratic portion of the measured ^gnals . 

Generally, a correlation canceler has two signal' inputs 
and one output. One of the inputs is either the secondary 
reference n' (t) or the primary reference s' (t) which are 
correlated, respectively, to the secondary signal portions 
n(t) and the primary signal portions s(t) present in the 
composite signal S(t) . The other input is for the composite 
signal S{t) . Ideally, the output of the corre-lsrticn cancs-l-cr 
s"(t) or n"(t) corresponds, respectively, -to -th^ , r^im?r y 
signal s(t) or the secondary signal n(t) portions only. 
Often, the most difficult task in the application of 
correlation cancelers is determining the reference sig^^ial-s- 
n' (t) and s' (t) which are correlated to the secondary n{t) 
and primary s(t) portions, respectively, of the measured 
signal S(t) since, as discussed above, these portions are 
quite difficult to isolate from the measured signal S(t) . In 
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the signal processor of the present invention, either a 
secondary reference n' (t) or a primary reference s' (t) is 
determined from two composite signals measured 
simultaneously, or nearly simultaneously, at two different 
5 wavelengths, Xa and Xb. 

A block diagram of a generic monitor incorporating a 
signal processor according to the present invention, and a 
correlation cancelisr- shown in Fiyuxes "4ci and— 4-b- 'Two^ 
measured signals, S;,^(t) and S;^(t), are acquired by a 

10 detector 20. One skilled in the art will realize that far 

some physiological measurements ^ more than cn«-, detector msf 
be advanta geous . Each signal is conditioned by a signal 
conditioner 22a and 22b. Conditioning includes, but is net: 
limited to, such procedures as filtering the signals to 

15 remove constant portions and amplifying the signals for ea se - 

of manipulation. The signals are then converted to digital 
data by an analog- to-digital converter 24a and 24b. The 
first measured signal S^^it) comprises a first primary signal 
portion, labeled herein S;^,(t), and a first secondary signal 

20 portion, labeled herein n^^Ct) . The second measured signal 

Sxb(t:) is at least partially correlated to the first measured 
signal S^^it) and comprises a second primary s-ignal portion-, 
labeled herein s^^^it) , and a second secondary signal portion, 
labeled herein n;aj(t)- Typically the first and second 

25 secondary signal portions, n^^Ct) and n;,b{t) , are uncorrelated 

and/or erratic with respect to the primary signal portions 
S3,^(t) and Sxb(t)' The secondary signal portions nj^^(t) and 
nxb{t) are often caused by motion of ,a patient in 
physiological measurements . 

30 The signals S^^{ti) and S^^i'c) are input to a reference 

processor 26. The reference processor 26 multiplies the 
second measured signal S.^Ct) by either a factor 
ra = s^^{t) /s^^it) or a factor r^ = n^^{M) /n^^it) and then- 
subtracts the second measured signal Sxb(t) from the first 

35 measured signal S^^it) . The signal coefficient factors r^ and 

r are determined to cause either the primary signal portions 
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S;^, (t) and S;^ij(t) or the secondary signal portions n^aCt) and 
nxb(t) to cancel, respectively, when the two signals S^^Ct) 
and Sxb(t) are subtracted. Thus, the output of the reference 
processor 26 is either a secondary reference signal n' (t) = 
5 n;,a(t:) - r^n^bCt), in Figure 4a, which is correlated to both of 

the secondary signal portions nx»(t) and nxb(t) or a primary 
reference signal s' (t) = Sx^Ct) - r^s^bCt) , in Figure 4b, which 
is correlated to both of the primary signal portions Sx^(t) 
and Sxb(t) . A reference signal n' (t) or s' it) d^-in?ut^ ^ "^rtr g 

10 with one of the measured signals Sx»(t) or S^bCt)/ to -=3- 

correlation canceler 27 which uses the reference signal n' (t) 
or s' (t) to remove either the secondary signal portions n^aCt) 
or nxb(t) or the primary signal portions s^aCt) or s^bCt) 
from the measured s-ign-a-l S^a(t) or Sxb(t} , Ths-cutput of the- 

15 correlation canceler 27 is a good primary signal 

approximation s" (t) or secondary signal approximation n" (t) . 
In one embodiment, the approximation s"(t) or n" (t) is 
displayed on a display 28. 

In one embodiment, an adaptive noise canceler 30, an 

20 example of which is shown in block diagram form in Figure. 5a, 

is employed as the correlation canceler 27, to remove either 
one of the erratic, secondary signal portions nxa(t) and 
n^bCt) from the first and second signals Sxa(t:) and Sxb(t:) . 
The adaptive noise canceler 30 in Figure 5a has as one input 

25 a sample of the secondary reference n' (t). which is correlated 

to the secondary signal portions nxa(t) and nxb(t) . The 
secondary reference n' (t) is determined from the two measured 
signals SxaCt) and Sxb(t) by the processor 26 ^the- pr-eeent 
invention as described herein. A second input to tixB 

30 adaptive noise canceler, is a sample of either the first -or 

second composite measured signals Sxa(t) = Sxa(t) + nxa(t) or 
SxbCt) = Sxb(t) + nxb(t) . 

The adaptive noise canceler 30, in Figure 5b, may also 
be employed to remove "either one of primary signal portions" 

35 Sxa{t) and s^^it) from the first and second measured signals 

Sxa(t) and Sxb(t) . The adaptive noise canceler 30 has as one 
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input a sample of the primary reference s' (t) which is 
correlated to the primary signal portions s^^Ct) and Sxb(t) . 
The primary reference s' (t) is determined from the two 
measured signals Sx.{t) and> Sxb(t) by the processor 26 of the 
5 present invention as described herein. A second input to the 

adaptive noise canceler 30 is a sample of either the first or 
second measured signals Sx.(t) = Sx.(t) + nx.(t) or Sx^dr) = 

Sxb(t) + nxb(t) . 

The adaptive noise canceler 30 functions to remove 
10 frequencies common to both the reference n' (t) or s' (t) and_ 

the measured signcrl Sx.(t) or Sxi,{t) . Since the reference- 
signals are correlated to either the secondary signal 
portions nx.(t) and nxb(t) or the primary signal portions- 
Sx,(t) and Sx^CU-.- the reference signals wilX. -bP... 
15 correspondingly erratic or well behaved. The adaptive noise- 

canceler 30 acts in a manner which may be analogized to a 
dynamic multiple notch filter based on the spectral 
distribution of the reference signal n' (t) or s' (t) . 

Figure 5c illustrates an exemplary transfer function of 
20 a multiple notch filter. The notches, or dips in. the 

amplitude of the transfer function, indicate frequencies 
which are attenuated or removed when a signal passe^s-throngh- 
the notch filter. The output of the notch filter is the 
composite signal having frequencies at which a notch is 
25 present removed. In the analogy to an adaptive noise 
canceler 30, the frequencies at which notches are present 
change continuously based upon the inputs to the adaptive 
noise canceler 30. • , , 

The adaptive noise canceler 30 (Figures 5 a and 5b) 
30 produces an output signal, labeled herein as s'\^{t) , 

s"xb(t), n"x.(t) or n"xb(t) which is fed back to an internal 
processor 32 within the adaptive noise canceler 30. The 
internal processor 32 automatically adjusts its own transfer 
function according to a predetermined algorithm such that the 
3 5 output of the internal processor 32 labeled bx(t) in Figure 

5a and Cx{t) in Figure 5b, closely resembles either the 
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secondary signal portion nx^Ct) or n^bCt:) or the primary 
signal portion s^^(t) or Sxb(t). The output h^{t) of the 
internal processor 32 in Figure 5a is subtracted from the 
measured signal, Sx*(t) or SxbCt) , yielding a signal output 
5 s"x^(t) = Sx*(t) + nx^(t) - bx^{t) or a signal output s"xb(t) = 

Sx5(t) +nxb(t) -bxb(t) . The internal processor optimizes s"xa(t) 
or s"xb(t:) such that s"x«(t:) or s"xb(t) is approximately equal 
to the primary signal s^*(t) or Sxb(t), respecti\':ely . The 
output Cx{t) of the internal processor 32 in Figure 5b xs* 

10 subtracted from the measured signal, S^^it) or Sxb(t), 

yielding a signal onfpnr given by n"x«{t) = Sx,(t) + - 
Cx^{t) or a signal outpu t given by n"xb(t) = Sxb(t) + nxbCt) - 
Cxb(t), The internal processor optimizes n"x^{t) or n"xb(t) 
such that n"x,(t) or n"xb(t) is approximately equal to the 

15 secondary signal portion nx^(t) or nxb(t) , respectively. 

One algorithm which may be used for the adjustment of 
the transfer function of the internal processor 32 is a 
least -squares algorithm, as described in Chapter 6 and 
Chapter 12 of the book Adaptive Signal Processing by Bernard 

20 Widrow and Samuel Stearns, published by Prentice Hall, 

copyright 1985. This entire book, including Chapters 6 and 
12, is hereby incorporated herein by reference. 

Adaptive processors 30 in Figures 5a and 5b have been 
successfully applied to a number of problems including 

25 antenna sidelobe canceling, pattern recognition, the 

elimination of periodic interference in general, and the 
elimination of echoes on long distance telephone transmission 
lines. However, considerable ingenuity is often required to 
find a suitable reference signal n' (t) or s' (t) since the 

30 portions nxa(t), nxb{t), Sxa(t) and Sxb(t) cannot easily be 

separated from the measured composite signals Sxa(t) and 
Sxb(t) . If either the actual secondary portion nxa(t) or 
nxb{t) or the primary signal portion Sxa(t) or Sxb(t) were a 
priori available, techniques such as correlation cancellation 

3 5 would not be necessary. 
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GENERALIZED DETERMINATION OF PRIMARY AND SECONDARY 

REFERENCE SIGNALS 
An explanation which describes how the reference signals 
n' (t) and s' (t) may be determined follows* A first signal is 
5 measured at, for example, a wavelength Xa, by a detector 

yielding a signal S^^Ct) : 

S^*(t) = Sx*(t) + nx^(t) (1; 

10 where Sx*(t) is the primary signal portion and Uj^Ct) is th^ 

secondary signal portion, 

A similar measurement is taken simultaneously, or nearly 
simultaneously, at a different wavelength, Xb, yielding: 

15 Sxfe(U = -S;,b(t) + n..t,{U. (2) 

Note that as long as the measurements, S;^^(t) and Sxb(t), are 
taken substantially simultaneously, the secondary signal 
components, n^.Ct) and nxb(t:), are correlated because any 
20 random or erratic functions affect each measurement in nearly 

the same fashion. The substantially predictable primary 
signal components, s^^^it) and Sxb(t) , are also correlated to 
one another. 

To obtain the reference signals n' (t) and s' (t) , the 
25 measured signals S^^(t) and Sxb(t) are transformed to 

eliminate, respectively, the primary or secondary signal 
components. In accordance with the present invention one wa^^ 
of doing this is to find proportionality constants, r^ and r^, 
between the primary signal portions Sxa(t) and Sxb(t) and che 
30 secondary signal portions nxa(t) and n^^it) such tliat ^tite * 

signals can be modeled as follows: 

Sxa(t) = r^ Sxb(t:) 

(3) 

n^Jt) = r, n,b{t) . 
3 5 In accordance with the inventive signal model of the present 

invention, these proportionality relationships can be 
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satisfied in many measurements, including but not limited to 
absorption measurements and physiological measurements. 
Additionally, in accordance with the signal model of the 
present invention, in most measurements, the proportionality 
5 constants r^ and r^ can be determined such that: 

nxa(t) 1^ r^n;,j,(t) 

(4) 

S;,.(t) ^ r^Sxb(t) • 
Multiplying equation {2) by r^ and then subtracting 
10 equation (2) from equation (1) results in a single equation- 

wherein the primary signal terms Sx»(t) and Sxb(t) ca n c e l r 

n' (t) = SxaCt) - r^S^i^it) = nx^(t) - r^n;,b(t) ; {5a} 

15 a non-zero g-i gna1 which is correlated to each secondary 

signal portion n;^a(t) and n^bC^) crHn used cis~l:iitr^ 

secondary reference n' (t) in a correlation canceler such as 
an adaptive noise canceler* 

Multiplying equation (2) by r^ and then subtracting 

20 equation (2) from equation (1) results in a single equation 

wherein the secondary signal terms nxa(t) aTid Tixb(t) cancel, 
leaving: 

s'{t) = Sxa(t) - S,b(t) - Sx,{t) - r, Sxb(t); (5b) 

25 

a non-zero signal which is correlated to each of the primary 
signal portions Sxa(t:) and Sxb(t) and can be used as the 
signal reference s' (t) in a correlation canceler such as an 
adaptive noise canceler. 

3 0 EXAMPLE OF DETERMINATION OF PRIMARY AND SECONDART 

REFERENCE SIGNALS IN AN ABSORPTIVE SYSTEM 
Correlation canceling is particularly useful in a large 
number of measurements generally described as absorption 
measurements. An example of an absorption type monitor wh i rrh- 

3 5 can advantageously employ correlation canceling, such as 

adaptive noise canceling, based upon a reference n' (t) or 
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s' (t) determined by a processor of the present invention is 
one which determines the concentration of an energy absorbing 
constituent within an absorbing material when the material is 
subject to change. Such changes can be caused by forces 
5 about which information is desired or primary, or 

alternatively, by random or erratic secondary forces such as 
a mechanical force on the material. Random or erratic 
interference, such as motion, generates secondary components 
in the measured signal > These secondary components can he 
10 removed or derived by the correlation canceler if a suitable^ 

secondary reFerence n' Ct) or primary reference s' (t) is 
known. 

A schematic N constituent absorbing material cor upxiblii^ 
a container having. H different abs-o>">^'' ^/-^nc^f--; ^MOTn^c:_„.. 

15 labeled A^, A2, A3,... A^, is shown in Figure 6'a. "Th-e 

constituents A, through A*, in Figure 6a are arranged in a 
generally orderly, layered fashion within the container 42. 
An example of a particular type of absorptive system is one 
in which light energy passes through the container 42 and is 

20 absorbed according to the generalized Beer-Lambert Law .of 

light absorption. For light of wavelength Xa, this 
attenuation may be approximated by: 



i-l 



25 Initially transforming the signal by taking -the nacural 

logarithm of both sides and manipulating terms, the signal io" 
transformed such that the signal components are combined hy 
addition rather than multiplication, i.e.: 



5,,=ln{Vl)=X;€,.,3C,x, (7) 
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where Iq is the incident light energy intensity; I is the 
transmitted light energy intensity; e^^xa is the absorption 
coefficient of the i*"^ constituent at the wavelength Xa; x^Ct) 
is the optical path length of i*"^ layer, i.e., the thickness 
5 of material of the i^^ layer through which optical energy 

passes; and c^Ct) is the concentration of the i^^ constituent 
in the volume associated with the thickness x^Ct)* The 
absorption coefficients through are known values which 
are constant at each wavelength. Most co ncern t rati . or r^{t) 

10 through c^Ct) are typically unknown, as axe most of trhe 

optical path lengths x^Ct) of each layer. The total optical 
path length is the sum of each of the individual optical path 
lengths Xi(t) of each layer. 

When the materxal is not subj-ect to any f crc-es— A^hz.ch 

15 cause change in the thicknesses of the layers, the optical 

path, length. jq:L .each—layer^. ^;-(.tX/ is generally constant. 
This results in generally constant attenuation of the optical 
energy and thus, a generally constant offset in the measured 
signal. Typically, this offset portion of the signal is of 

20 little interest since knowledge about a force which perturbs 

the material is usually desired. Any signal portion outside 
of a known bandwidth of interest, including the constant 
undesired signal portion resulting from the generally 
constant absorption of the constituents when not subject to 

25 change, is removed. This is easily accomplished by 

traditional band pass filtering techniques. However, when 
the material is subject to forces, each layer of constituents 
may be affected by the perturbation differently than other 
layers. Some perturbations of the opt i c al path len gt h s err 

30 each layer Xi(t) may result in excursions in the measured 

signal which represent desired or primary information. Otheir 
perturbations of the optical path length of each layer Xi-(t) 
cause - undesired or secondary excursions which mask primary 
information in the measured signal • Secondary sigrrarl 

35 components associated with secondary excursions must also be 

removed to obtain primary information from the measured 
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signal. Similarly, the ability to compute secondary signal 
components caused by secondary excursions directly allows one 
to obtain primary signal components from the measured signal 
via simple subtraction, or correlation cancellation 
5 techniques . 

The correlation canceler may selectively remove from the 
composite signal^ measured after being transmitted through or 
reflected from the absorbing material, either- the spcnndary^ 
or the primary signal components caused by forces 

10 perturb or change the material differently from the forces 

which perturbed or changed the material to cause- 
respectively, either the primary or secondary signal 
component. For the purposes of illustration, it will he- 
assumed that the portion of the measured signal which is 

15 deemed to be the primary signal s^^Ct) is the attenuation term 

esCgXsCt) associated with a constituent of interest, namely 
A5, and that the layer of constituent A5 is affected by 
perturbations different than each of the layers of other 
constituents through and A^ through Ajj. An example of 

20 such a situation is when layer A5 is subject to forces about 

which information is deemed to be primary and, additionally/ 
the entire material is subject to forces which affect each of 
the layers. In this case, since the total force affecting 
the layer of constituent A5 is different than the total 

25 forces affecting each of the other layers and information is 

deemed to be primary about the forces and resultant 
perturbation of the layer of constituent Ag, attenuation 
terms due to constituents A^ through A^ and Ag through make 
up the secondary signal portion n^^it) . 'Even if the 

30 additional forces which affect the entire material cause the 

same perturbation in each layer, including the layer of Ag, 
the total forces on the layer of constituent Ag cause it to 
have different total perturbation than each of the other 
layers of constituents A^ through A4 and A^ through A^^,, 

35 It is often the case that the total perturbation 

affecting the layers associated with the secondary signal 
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components is caused by random or erratic forces. This 
causes the thickness of layers to change erratically and the 
optical path length of each layer, x^Ct), to change 
erratically, thereby producing a random or erratic secondary 
5 signal component n^^Ct) • However, regardless of whether or 

not the secondary signal portion nj^^(t) is erratic, the 
secondary signal component nx, (t) can be either removed or 
derived via a correlation canceler, such as an adaptive noise 
canceler, having as one input, respectively, a secondaxy* 

10 reference n' (t) or a primary reference s' (t) determined by a 

processor of the present invention as long as the 
perturbation on layers other than the layer of constituent Ag 
is different than the perturbation on the layer of 
constituent Ag. "The correlacion canceler yields a good 

15 approximation to either the primary signal Sxa(t) or the 

to the primary signal is obtained, the concentration of the 
constituent of interest, CjCt), can often be determined since 
in some physiological measurements, the thickness of the 

20 primary signal component, X5(t) in this example, is known or 

can be determined. 

The correlation canceler utilizes either the secondary 
reference n' (t) or the primary reference s' (t) determined 
from' two substantially simultaneously measured signals Sxa(t) 

25 and Sxb(t), Sx^^Ct:) is determined as above in equation (7). 

Sxb(t) is determined similarly at a different wavelength Xb, 
To find either the secondary reference n' (t) or the primary 
reference s' (t) , attenuated transmitted energy l*s -measured ac 
the two different wavelengths Xa and Xb and transformed via 

30 logarithmic conversion. The signals Sxa(t) and Sx^Ct) can 

then be written (logarithm converted) as: 

4 S 



Sxa 0=€5,,,C5Xs{t)+nxa(t) 
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(9) 



4 ^ 
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Sxb(t)=C5.XbCs3C5(t)+nxb(t) (11) 

Further transformations of the signals are the 
5_ proportionality relationships in accordance with the signal 
model of the present invention defining r, and r. . similar tn . 
equation (3) , which allows determination of a nx3ise--ref-er-eTi,c£- 
n' (t) and a primary reference s' (t) . These are-: 

€s.x. = r,6s.xb ^12a) 
10 nx. = r^xb <l2b) 

where 

€s. xa ^ €s, xb- ^13b) 

j;^3 Otten uHS ^^cxSo w^-/v-**-^>iWc;N,-..w**w X » *■ 

15 simultaneously satisfied. Multiplying equation (11) by r, 
and subtracting the result from equation (9) yields a non- 
zero secondary reference which is a linear sum of secondary 
signal components: 

n' (t) = S^^[t) - r.Sxb(t) = nxa(t) - r,n,^(t) (14a) 



20 



i 6i,x.CiXi(t) + I ei.x.CiXi(t) 



Z r^ei,xbCiXi{t)+ i r,ei,xbCiXi(t) 



(15a) 



= i CiXi(t) t Ci.xa - r^ei.xb 3 + 



25 



i CiXi(t) ( €i.x4 - r,€i,^i, ] 



(16a) 
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Multiplying equation (11) by and subtracting the 
result from equation (9) yields a primary reference which is 
a linear sum of primary signal components: 



10 



s' (t) - S;,^(t) - r^xb(t:) = SA*(t) - r,sxb(t) (14b) 

= CsXs(t)€s,x* - C5X5(t)€s,xb (15b) 

= CsXs(t) [€s,x* " €s,xbj. -U^- 

A sample of either the secondary reference n' (t) or the 
primary reference s' (t) , and a sample of either measured 
signal Sxa(t) or Sxb(t), are input to a correlation canceier 



27/ such as "arr* adapt iVo noise cance-l-^r 30/ ejcsntplc cf 
15 which is shown in Figures 5a and 5b and a preferred example 

of which is /-^^ o^^i<=gg>H h^yi=>Tn landfir th.e .heading PREFERRED 
CORRELATION CANCELER USING A JOINT PROCESS ESTIMATOR 
IMPLEMENTATION • The correlation canceier 27 removes either 
the secondary portion nx*(t) or nxb(t), or the primary 
20 portions, s^^it) or Sxb(t), of the measured signal yielding a 

good approximation to either the primary signals s'\^.{t) - 
6s^xaCsX5(t) or s"xi,(t:) - €s,xbCsX5(t) or the secondary signals 
n"xa(t) - nxa(t) or n"xb(t) - nxb(t) . In the event that the 
primairy signals are obtained, the concentration Cs{t) may 
25 then be determined from the approximation to the primary 

signal s«x^(t) or s"xb(t) according to: 

C5(t) - s"x,(t)/€s,xaX5(t:) , (17a) ^ 

or 

30 Cs(t) « s"xb{t)/6s,xbXs(t) (12d) 

As discussed previously, the absorption coefficients are 
constant at each wavelength Xa and Xb and the thickness— of- 
the primary signal component, Xs(t) in this example, is often 
known or can be determined as a function of time, thereby 

35 allowing calculation of the concentration c^{t) of 

constituent A5. 
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DETERMINATION OF CONCENTRATION OR SATURATION 



IN A VOLUME CONTAINING MORE THAN ONE CONSTITUENT 
Referring to Figure 6b, another material having N 
different constituents arranged in layers is shown , In this 
5 material/ two constituents Ag and A^ are found within one 

layer having thickness X5^fi{t) = XjCt) + XgCt), located 
generally randomly within the layer. This is analogous to 
combining the layers of constituents 7^ and^A^ in Figure ea, 
A combination of layers , such as the combination of layers of 
10 constituents Ag and A^, is feasible when the two layers ^r-*^ 

under the same total forces which result in the same change 
of the optical pach lengths XjCt) and x^Ct) of che layers. 

Often it is desirable to find the concentration or LLc 

saturation^ -j ^ - ^ p^yrf^nt- nnnn^nt-y^iti nn ^ nf on^^ r-on gf- -i Anh 

15 within a given thickness which contains more than oner 

constituent and is subject to unique forces, A determination 
of the concentration or the saturation of a constituent 
within a given volume may be made with any number of 
constituents in the volume subject to the same total forces 

20 and therefore under the same perturbation or change. To 

determine the saturation of one constituent in a volume 
comprising many cuay L i L uent s , as many measuxt:d' slgixcils as~" 
there are constituents which absorb incident light energy are 
necessary. It will be understood that constituents which do 

25 not absorb light energy are not consequential in the 

determination of saturation. To determine the concentration, 
as many signals as there are constituents which absorb 
incident ligkt energy are necessary as weUL a-s i nf or-m^fi on 
about the sum of concentrations. 

3 0 It is often the case that a thickness under unique 

motion contains only two constituents. For example, it may 
be desirable to know the concentration or saturation of Ac 
within a given volume which contains A^ and A^ . In this case, 
the primary signals S;^3(t) and Sxb{t) comprise terms related 

3 5 to both As and A^ so that a determination of the concentration 

or saturation of As or A^ in the volume may be made. A 
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determination of saturation is discussed herein, it will be 
understood that the concentration of A5 in a volume 
containing both and could also be determined if it is 
known that Ag+A^sl, i.e., that there are no constituents in 
5 the volume which do not absorb incident light energy at the 

particular measurement wavelengths chosen. The measured 
signals Sxa(t) and S^bCt) can be written (logarithm converted) 
as : 

10 Sy^it) = €5,xaCsX5.6{t) + ^ ^^^^c^K^^^ {t) + xi^it) CXFa) 

= Sxa(t) + nxa(t) (18b) 



15 



It is also often the case that there may be two or more 
thicknesses within a medium each containing the same two 

20 constituents but each experiencing a separate motion as -in 

Figure 6c. For example, it may be desirable to know the 
concentration or saturation of Ag within a given volume which 
contains Ag and A^ as well as the concentration or saturation 
of A3 within a given volume which contains A3 and A^, A3 and 

25 A4 having the same constituency as A5 and A^ respectively. In 

this case, the primary signals Sxa(t) and Sxb(t) agaiir 
comprise terms related to both Ag and Ag and portions of the 
secondary signals nxa(t) and n;^bt^) comprise cerms relaced zo 
both A3 and A^ . The layers, A3 and A^, do not enter into trher 

3 0 primary equation because they are assumed ta fae perturbed'-^by 

a different frequency, or random or erratic secondary forces 
which are uncorrelated with the primary force. Since 
constituents 3 and 5 as well as constituents 4 and 6 are 
taken to be the same, they have the same absorption 

35 coefficients (i.e., Cj^^a = ^g^x*; ^i^xh = ^s^xb; x* = ^g^xa 

^4,xb - ^6 xb- Generally speaking, however, A3 and A^ will have 
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different concentrations than Ag and and will therefore 
have a different saturation. Consequently a single 
constituent within a medium may have one or more satxirations 
associated with it. The primary and secondary signals 
5 according to this model may be written as: 

Sx*(t) = te5,x*Cs+ €fi,x*Cfi ] Xs,s(t) {20a) 
10 + X fi,x*CiXi(t)+ f ei,x*CiXi(t) (20b) 

i«l i«7 



15 



2 * 

+ 1 €i,xb CiXi{t)+ X ei,xbCiXi(t) • (21b) 



20 nxb(t) = tes^xbCj + e«,xbC4] X3,4(t) + n^{t) (21c) 

where signals nxa(t) and nxb(t) are similar to the secondary 
signals nxa(t) and nxb(t) except for the omission of the 3, 4 
layer. ^ ' - 

25 Any signal portions whether primairy or secondary, 

outside of a known bandwidth of interest, including the 
constant undesired secondary signal portion resulting from 
the generally constant absorption of the constituents when 
not under perturbation, should be removed to determine an 

30 approximation to either the primary signal or the secondary 

signal within the bandwidth of interest. This is easily 
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accomplished by traditional band pass filtering techniques. 
As in the previous example, it is often the case that the 
total perturbation or change affecting the layers associated 
with the secondary signal components is caused by random or 
5 erratic forces, causing the thickness of each layer, or the 

optical path length of each layer, x^Ct), to change 
erratically, producing a random or erratic secondary signal 
component nx^Ct) . Regardless of whether or not the secondary/ 
signal portion n^^Ct) is erratic, the secondary signal 

10 component n^aCt) can be removed or derived via a correlation 
canceler, such as an adaptive noise canceler, having as one 
input a secondary reference n' (t) or a primary reference 
s' (t) determined by a processor of the present invention as 
long as the perturbation in layers other than tne layer or 

15 constituents A5 and is different than the perturbation in 

the layer of 'coijcstituernts A^ and A^ . 2 jL the^r t h^ Sjrr a t c 
secondary signal components n;^a(t) and n^^^it) or the primary 
components Sx^Ct) and s^^{t) may advantageously be removed 
from equations (18) and (19), or alternatively equations (20) 

20 and (21) , by a correlation canceler. The correlation 

canceler, again, requires a sample of either the primary 
reference s' (t) or the secondary reference n' (t) and a sample 
of either of the composite signals Sxa(t) or Sxb(t) of 
equations (18) and (19) • 

25 DETERMINATION OF PRIMARY AND SECONDARY REFERENCE SIGNALS 

FOR SATURATION MEASUREMENTS 
One method for determining reference signals s' (t) or 
n' (t) from the measured signals Sxa{t) and S^bCt) in 
accordance with one aspect of the invention is what will be 

30 referred to as the constant saturation approach. In this 

approach, it is assumed that the saturation of A5 in the 
volume containing A5 and A^ and the saturation of A3 in the 
volume containing A3 and A4 remains relatively constant over 
some period of time, i:e,: 

35 

Saturation (A5 (t) ) = Cs(t) / [Cs{t) + c^(t)] (22a) 
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Saturation {A3 (t) ) = c^it) / [c^Ct) + c,(t)] {22b) 

Saturation (As {t) ) = {l + [c^ {t) /Cj {t) ] {23a) 

5 Saturation (A3 ( t) ) = {l + Cc^ (t) /c^ (t) ] (23b) 

are substantially constant over many samples of the measured 
signals Sx^ and S^^. This assumption is accurate over many 
samples since saturation generally changes rela^tively slo-.ciy 
10 in physiological systems - 

The constant saturation assumption is equivalent to 
assuming that: 



15 



30 



r-.jf'\J,r- Jf-^ = const-ant— i2Ah)^ 



since the only other term in equations (23a) and (23b) is a 
constant, namely the numeral 1. 
20 Using this assumption, the proportionality constants, r, 

and r^ which allow determination of the secondary reference 
signal n' (t) and the primary reference signal s' (t) in the 
constant saturation method are: 



25 Cs.xa Cs Xs,s(t) + 6«,x. Cs Xs.s(t) 



(25a) 



^S.Xb Xs,j(t) + €«,xb X5,s(t) 

= s,,{t)/s,,{t) • (26a) 



(27a) 



35 Cs.xa (Cs/c«) + €«.x. 



(28a) 
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s"x3i (t) /s*\b(t) = constantj; where (29a} 
nx*(t:) ^ r^(t)n;,b(t:) (30a} 

and 



10 



20 



= nxJt)/n;^Ct) (26b) 



15 €s,xa(C3/C4) + 



Xa 



(28b) 



- n"xa(t) /^"xb(t^) = constant^; where (29b) 



In accordance with the present invention, it is often 
the case that both equations (26) and (30) can be 

25 simultaneously satisfied to determine the proportionality 

constants r^ and r^. Additionally, the absorption 

coefficients at each wavelength e^^y^, e^^, and e^^xb ^^re 

constant and the central assumption of 'the constant 
saturation method is that Cs(t)/c6(t) and C3(t}/c4(t) are 

30 constant over many sample periods. Thus, new proportionality 

constants r^ and r^ may be determined every few samples from 
new approximations to either the primary or secondary signal 
as output from the correlation canceler. Thus, the 
approximations "to either the primary signals Sxa(t) and Sxb(^) 

35 or the secondary signals nxa(t) and nxb(t) , found by the 

correlation canceler for a substantially immediately 
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preceding set of samples of the measured signals Sx^(t) and 
SxbCt) are used in a processor of the present invention for 
calculating t^he proportionality constants, r^ and r^, for the 
next set of samples of the measured signals Sx^Ct) and Sxb(t) • 
5 Multiplying equation (19) by and subtracting the 

resulting equation from equation (18) yields a non-zero, 
secondary reference signal: 

n' (t)=Sx*(t)-r.S,b{t)-n,^(t)-r^nx5(t) • (31a) 

10 

Multiplying equation (19) by r^ and subtracting the 
resulting equation from equation (18) yields a non-zero 
primary reference signal: 

15 ^' (t) -S.,(t) -x^xbCt:) =Sxa(t:) -r^xi.(^^-- , 

When using the constant saturation method in patient 
monitoring, initial proportionality coefficients can be 
determined as further explained below. It is not necessary 
20 for the patient to remain motionless even for- an 

initialization period. With values for the proportionality 
coefficients r^ and r^ determined, a correlation canceler -may- 
be utilized with a secondary reference n' (t) or a primary 
reference s' (t) . 
25 DETERMINATION OF SIGNAL COEFFICIENTS FOR 

PRIMARY AND SECONDARY REFERENCE SIGNALS 
TT.qiNG THE CONSTANT SATURATION METHOD 
In accordance with one aspect of the present invention, 
the reference processor 26 of Figures 4a and Figure 4b of the 
3 0 present invention may be configured to multiply the second- 

measured assumed signal Sxb(t:) =Sxb(t) +n;^b(t) by each of a 
plurality of signal coefficients r^, rj, • - • r^, and then 
subtract each result from the first measured signal 
Sxa{t:) =Sxa(t:) +nx. (t) to obtain a plurality of reference signals- 

35 

R' (r, t) = Sxa(t) - r sxb(t) + n^Jt) - r nx^Ct) (32) 
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for r = r^, • • • as shown in Figure 7a. In other words, 

a plurality of signal coefficients are chosen to represent a 
cross section of possible signal coefficients. 

In order to determine either the primary reference s' (t) 
5 or the secondary reference n' (t) from the above plurality of 

reference signals of equation (32) , signal coefficients r^ 
and r^ are determined from the plurality of assumed si gn a l 
coefficients r^, r^, — • r„. The coef ficienJts jr. :an± r , -rn^ 
selected such that they cause either the primary signal 
10 portions s^^Ct) and s^bCt) or the secondary signal p ort inns 

n^.Ct) and n;^b{t) to cancel or nearly canceL wlien tJiey are 
substituted into the reference function RMr, t) , e, g, 

Sx*(t) = r.Sxi>(t) {33a) 

15 

n' (t) = R' (r., t) = nxa(t) - r^nxb(t) (33c) 

20 s' (t) = R' (r^, t) = s,,(t) - r,s,b(t). (33d) 

In other words, coefficients r^ and r^ are selected at 
values which reflect the minimum of correlation between the 
primary signal portions and the secondary signal portions. - 

25 In practice, one does not usually have significant prior 

information about either the primary signal portions S;^^(t) 
and Sxb(t) or the secondary signal portions nxa(t) and n^t^^) 
of the measured signals S^^it) and Sxb(t) . The -lack of this 
information makes it difficult to determine which of rhe 

30 plurality of coefficients r^, r^ correspond to the^ 

signal coefficients r^ = s^Jt) /s^^it) and r^ = n^^it) /n^^it.) . 

One approach to determine the signal coefficients r, and" 
■r^ from the plurality of coefficients r^, rj, — • reemploys 
the use of a correlation canceler 27, such as an adaptive 

35 noise canceler, which takes a first input which corresponds 

to one of the measured signals S;^, (t) or S^^^it) and takes a 
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second input which corresponds to successively each one of 
the plurality of reference signals R' (r^, t) , R' (rj, t) , — 
R' (r^,, t) as shown in Figure 7a, For each of the reference 
signals R' (r^, t) , R' (rjr t) , • • • , R' (r^,, t) the corresponding 
5 output of the correlation canceler 27 is input to a "squares" 

operation 28 which squares the output of the correlation 
canceler 27, The output of the squares operation 28 ±s_ 
provided to an integrator 29 for forming a cumulative output 
signal (a summation of the squares) . The cumulative o utput • 

10 signal is subsequently input to an extremum detector 31. The 

purpose of the extremum detector 31 ^ is to chose signal 
coefficients r^ and r^ from the set r^, x^, • • • i^n t>y observing 
which provide a maximum in the cumulative output signal as in 
Figures 7b and To. Tn other words, coefficients wnich 

15 provide a maximum integrated output, such as energy or power, 

from the correlation -canceler 27 corr-espond to the — signal- 
coefficients r^ and r^ which relate to a minimum correlation 
between the primary signal portions and the secondary signal 
portions in accordance with the signal model of the present 

20 invention. One could also configure a system geometry which 

would require one to locate the coefficients from the set r^, 
X2, * • • r^^ which provide a minimum or inflection in the 
cumulative output signal to identify the signal coefficients 
and r^, 

25 Use of a plurality of coefficients in the processor of 

the present invention in conjunction with a correlation, 
canceler 27 to determine the signal coefficients r^ and r^ may 
be demonstrated by using the properties of ^correlation 
cancellation. If x, y and 2 are taken to be any collection 

3 0 of three time varying signals, then the properties of some 

correlation cancelers C{x, y) may be defined as follows: 

Property (1) C(x, y) = 0 for y correlated (34a) 
Property (2) C{x, y) = x for x, -y uncorrelated (34b) 
35 Property (3) C(x + y, 2) = C(x, 2) + C{y, 2) {34c) 
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With properties (1), (2) and (3) it is easy, to demonstrate 
that the energy or power output of a correlation canceler 
with a first input which corresponds to one of the measured 
signals Sx*(t) or Sxb(t:) and a second input which corresponds 
5 to successively each one of a plurality of reference signals 

R' (r^, t), R' (rj, t) , R' (r^, t) can determine the signal 

coefficients r.- and r^ needed to produce the primary - r e f e r encB:: 
s' (t) and secondary reference n' (t) • If we take as a first 
input to the correlation canceler the measured signal S;^^(t) 
10 and as a second input the plurality of reference signals 

R' (ri, t), R' Cr^, U ^ R' {r„, t) then the outputs of the 

correlation canceler C{S;^(t), RMr^^t)) for j = 1, 2, •••^ n 
may be written as 

15 C(Sxa(t) + nx^(t),S;,^(t}- r^SxbCt) + nx,(t)- rjn;,b(t)) (35) 

where j = 1, 2, n and we have used the expressions 

R' (r, t) = S,Jt) - rSxb(t) (36) 

20 

SxaCt) = Sxa(t) + nx^(t) (37a) 

Sxb(t) = Sxb(t) + nxb(t), (37b) 

25 The use of property (3) allows one to expand equation (35) 

into two terms 

C(Sxa(t) .R' (r,t) )=C(Sxa(t:) .Sx,(t) -rSxb(t:) +nx^(t) -rnxb(t) ) 

30 +C(nx,(t) ,Sxa(t) -rsxb(t)+nx,(t)-rnxb{t)) (38) 

so that upon use of properties (1) and (2) the correlation 
canceler output is given by 

C(Sx,(t), R'"(rj,t))= Sxa(t)6(rj-rJ ^ n^, (t ) 6 (r^-rj (39) 

35 

where 6 (x) is the unit impulse function 
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6{x) = 0 if X 7i 0 
6{x) = 1 if X = 0. 
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5 The time variable, t, of the correlation canceler output 

CCSxaCt) , R' (r^, t)) may be eliminated by computing its energy 
or power. The energy of the correlation canceler output is 
given by 

10 Ex^(rj)=:/c^(Sx*(t) ,R' {rj,t)dt 

=5 (r^-r J fs\^{t) dt+5 (r^-rj Jn\^ (t ) dt . (41a) 

It should be understood that one could, equally well, 
have chosen the measured signal S^bCt:) ccs^ Llib f juj-st iTtptrt to 
15 the correlation canceler and the plurality of reference 

signals R' (r^, t) , R' {r^, t) , • - • . J?jL.(x^_U fhfi F^pnonrl , 

input- In this event, the correlation canceler energy output 
is 

20 Exb(r^)=/c^(Sxb{t) ,R' (r,t)dt Ulb) 

=6 (r^-rj /s\b (t) dt+6 (r^-r^) Sn\^{t) dt , 

It should also be understood that in practical 
situations the use of discrete time measurement signals may 

25 be employed as well as continuous time measurement signals, 

A system which performs a discrete transform (e.g., a 
saturation transform in the present example) in accordance 
with the present invention is described leirth , ir^f ererrce- to 
Figures 11-22, In the event that discrete time measureinerLt^ 

30 signals are used, integration approximation methods such as 

the trapezoid rule, midpoint rule. Tick's rule, Simpson's 
approximation or other techniques may be used to comput-e- the 
correlation canceler energy or power output. In the discrete 
time measurement signal case, the energy output of rhe 

3 5 correlation canceler may be written, using the trapezoid 

rule, as 
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Ex. 



(rj = 6(r.-r.)At{ f s\Jt,) - 0.5(s\.(to) + s^,(t„))} 



+6{r^-rjAt{ I n^^(ti)-0.5{n^.{to)+n^^(t„)) } (42a) 



i-O 



(r) = 6(r.-rjAt{ f s\b(ti) - 0,5{s\b(to) + s\i,{t„))} 



1«0 



+6(r.-rjAt{ I n\b(ti)-0.5{n\b(to)+n\b(tn)) } (42b) 



i-O 



where is the discrete time, tg is the initial time, t„ 
is the final time and At is the time between discrete time 

10 measurement samples. 

The energy functions given above, and shown in Figure 
7b, indicate that the correlation canceler output is usually 
zero due to correlation between the measured signal S^^Ct) or 
Sxb(t) and many of the plurality of reference signals 

15 R' (ri, t), R' (rj, t) , R' (r„, t) . However, the energy 

functions are non zero at values of r^ which correspond to 
cancellation of either the primary signal portions Sx.(t) and 
Sxb(t) or the secondary signal portions n^,. (t) and nxb(t) in 
the reference signal R' (r^, t) . These values ,co;rrespond to 

20 the signal coefficients r, and r^. 

It should be understood that there may be instances in 
time when either the primary signal portions S;,,(t) and s^b^^) 
or the secondary signal portions n^^{t) and n^bCt) are 
identically zero or nearly zero. In these cases, only one 

2 5 signal coefficient value will provide maximum energy or power 

-output of the correlation canceler. 
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Since there may be more than one signal coefficient 
value which provides maximum correlation canceler energy or 
power output, an ambiguity may arise. It may not be 
immediately obvious which signal coefficient together with 
5 the reference function R' (r, t) provides either the primary 

or secondary reference. In such cases, it is necessary to 
consider the constraints of the physical system at hand. For 
example, in pulse oximetry, it is known that arterial blood, 
whose signature is the primary plethysmographic wave, has. 

10 greater oxygen saturation than venous blood, whose signature- 

is the secondary erratic or random signal. Consequently, in 
pulse oximetry, the ratio of the primary signals due to 
arterial pulsation r. = S;,^ (t) /s^b (t) is the smaller of the two 
signal coefficient values while the ratio of the s^cor*<^?r^ ^ 

15 signals due to mainly venous blood dynamics r^ = n^a (t) /n^b (t) 

is the larger of the two signal coefficient values, assuming 
Xa = 660 nm and Xb = 910 nm. 

It- should also be understood that in practical 
implementations of the plurality of reference signals and 

20 cross correlator technique, the ideal features listed, as 

properties (1) , (2) and (3) above will not be precisely 
satisfied but will be approximations thereof. Tnerefore, ±rr 
practical implementations of this embodiment of the present 
invention, the correlation canceler energy curves depicted in 

25 Figure 7b will not consist of infinitely narrow delta 

functions but will have finite width associated with them as 
depicted in Figure 7c. 

It should also be understood that it is possible to ha^^ 
more than two signal coefficient values which produce maximum 

3 0 energy or power output from a correlation canceler. This. 

situation arises when the measured signals each contain more 
than two components each of which are related by a ratio as, 
follows : 
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S,,(t) = I fx..i{t) 

i*0 



(43) 



where 



Thus, reference signal techniques together with a 

10 correlation CcincellaJiion, such as an adaptive noise canceler. 

can be employed to decor-.pose a signal into two or niors signal 
components each of which is related by a ratio. 

PREFERRED CORRELATION CANCELER USING A 
■TOTNT PROCESS ESTIMATOR IMPLEMENTATION 

15 Once either the secondary reference n' (t) or the primary 

reference s' (t) is determined by the processor of the pr'esent 
invention, the correlation canceler can be implemented in 
either hardware or software. The preferred impl-ementation. of 
a correlation canceler is that of an adaptive noise canceler 

20 using a joint process estimator. 

The least mean squares (LMS) implementation of the 
internal processor 32 described above in conjunction with the 
adaptive noise canceler of Figures 5a and Figure 5b is 
relatively easy to implement, but lacks 'the speed at 

25 adaptation desirable for most physiological monitoring 

applications of the present invention. Thus, a fasrer- 
approach for adaptive noise canceling, called a least -squar=es- 
lattice joint process estimator model, is us^d. in one 
embodiment. ^ A joint process estimator 60 is shown 

30 diagrammatical ly in Figure 8 and is described in detail in 

Chapter 9 of Ada ptive Filter Theory by Simon Haykin, 
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published by Prentice-Hall, copyright 1986 • This entire 
book, including Chapter S, is hereby incorporated herein by 
reference • 

The function of the joint process estimator is to 
5 remove either the secondary signal portions n^^Ct) or n;^{t) 

or the primary signal portions s^aCt) or s^^it) from the 
measured signals Sy^(t) or S^^it) , yielding either a primary 
signal approximation s"xa(t) or s»^{t) or a secondary signal 
approximation n«x,{t) or n\^{t) . Thus, the joint prc c ™- 

10 estimator estimates either the value of the primary signsis^ 

S;,a(t) or Sxb(t) or the secondary signals n^^it) or n^bCt) . The 
inputs to the joint process estimator 60 are either the 
secondairy reference n' (t) or the primary reference s' (t) ana 
the compos it er-me€:sur^ed signal S^a't) cr S^b't) • Th^— cutpUw-X3 - 

15 a good approximation to the signal S^^it) or Sxb(t) with 

either the secondary signal or the primary signal removed, 
i.e. a good approximation to either s^^it) , s^^it) , n^aCt) 
or n^bCt) . 

The joint process estimator 60 of Figure 8 utilizes, in 

20 conjunction, a least square lattice predictor 70 and a 

regression filter 80. Either the secondary reference n' (t) 
or the primary reference s' (t) is input to the least square 
lattice predictor 70 while the measured signal Sxa(t) or 
Sxb(t) is input to the regression filter 80. For simplicity 

25 in the following description, Sy^it) will be the measured 

signal from which either the primary portion Sxa(t) or the 
secondary portion n^^^it) will be estimated by the joint 
process estimator 60, However, it will be r.ct-ed that Sxb(t) 
could also be input to the regression filtpr 8^ ^^^^ 

30 primary portion Sxb(t) or the secondary portion n^{t) of this 

signal could be estimated. 

The joint process estimator 60 removes all frequencies- 
that are present in both the reference n' (t) or s' (t) , and 
the measured signal S;,^(t). The secondary signal poxtrion- 

35 nxa(t) usually comprises frequencies unrelated to those of the 

primary signal portion S;,^(t) . It is improbable that the 
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secondary signal portion n^{t) would be of exactly the same 
spectral content as the primary signal portion S;,, (t). 
However, in the unlikely event that the spectral content of 
Sy^it) and nx.(t) are similar, this approach will not yield 
5 accurate results. Functionally, the joint process estimator 

60 compares the reference input signal n' (t) or s' (t) , which 
is correlated to either the secondary signal portion nx»(t) 
or the primary signal portion S;,.(t), and input signal S-^it) 
and removes all frequencies which are identic!. Trrcrer, xhB" 
10 joint process estimator 60 acts as a dynamic multiple notch 
filter to remove those frequencies in the secondary signal - 
component n;,,(t) as they change erratically with the motion 
of the patient or those frequencies in the primary signal 
component S;^(t) as they change with the arterial pulsacion 
15 of the patient. This yields a signal having substantially 
the same speuLxal cuiiLeut -3iTd-aTnpliti3t5e-a5-ei-t-h6-r-th€-pri~.ary - 
signal Sx,(t) or the secondary signal nx.(t) . Thus, the 
output s"x.(t) or n"x«(t) of the joint process estimator 60 is 
a very good approximation to either the primary signal Sxa(t) 
20 or the secondary signal nx,{t) . • • . 

The joint process estimator 60 can be divided xntc 
stages, beginning with a zero- stage and terminating in an m*^**- 
stage, as shown in Figure 8. Each stage, except for the 
zero-stage, is identical to every other stage. The zero- 
25 stage is an input stage for the joint process estimator 60. 

The first stage through the m'^^-stage work on the signa-1 
produced in the immediately previous stage, i.e., the (m-l)"''- 
stage, such that a good primary signal approximation s^xa^'t) 
or secondary signal approximation n"x,(t) is produced as 
30 output from the m^**- stage. 

The least-squares lattice predictor 70 comprises 
registers 90 and 92, summing elements 100 and 102, and delay 
elements 110. The registers 90 and 92 contain multiplicative, 
values of a forward reflection coefficient rf,„(t) and a- 
35 - backward reflection coefficient T^,„it) which multiply the 
reference signal n' (t) or s' (t) and signals derived from the 
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reference signal n' (t) or s' (t) . Each stage of the least- 
squares lattice predictor outputs a forward prediction error 
f„(t) and a backward prediction error b„(t) . The subscript 
m is indicative of the stage. 

For each set of samples, i.e. one sample of the 
reference signal n' (t) or s' (t) derived substantially 
simultaneously with one sample of the measured signal 
S^it) , the sample of the reference signal n' (t) or s' (t) is 
input to the least-squares lattice predictor 70, Zhe..z«rc-- 
stage forward prediction error fo(t) and the zero-st«ge- 
backward prediction error bo(t) are set equal to the 
reference signal n' (t) or s' (t> . The backward prediction 
error boCt) is delayed by one sample period by the delay 
element 110 in the first stage of the 1-ess^t— scjuar-ss l£tt;.c£- 
predictor 70. Thus, the immediately previous value of the 
xaf erence n' (ii) Qr_s.'_(t) is used in calculations involving 
the first-stage delay element 110. The zero-stage forward 
prediction error is added to the negative of the delayed 
zero-stage backward prediction error bo(t-l) multiplied by 
the forward reflection coefficient value rf,i(t) register. 90 
value, to produce a first-stage forward prediction error 
fi(t) . Additionally, the zero-stage forward prediction error 
fo(t) is multiplied by the backward reflection coefficient 
TbiCt) register 92 value and added to the delayed zero-stage 
backward prediction error bo{t-l) to produce a first-stage 
backward prediction error bi{t). In each subsequent stage, 
m, of the least square lattice predictor 70, the previous 
forward and backward prediction error values, f„.i(t) and b„. 
i(t-l), the backward prediction error be±ag delayed- iy isns. 
sample period, are used to produce values of the forward and 
backward prediction errors for the present stage, f„(t) and 
b„(t}. 

The backward prediction error b„(t) is fed to the 
concurrent stage, m, of the regress.ion filter 80. There it 
is input to a register 96, which contains a multiplicative 
regression coefficient value K„.x.(t) . For example, in the 



zero-stage of the regression filter 80, the zero-stage 
backward prediction error bo(t) is multiplied by the zero- 
stage regression coefficient <o,xa(t) register 96 value and 
subtracted from the measured value of the signal S^^Ct) at a 
5 summing element 106 to produce a first stage estimation error 

signal ei^x*(t)* The first-stage estimation error signal 
^i.x«(t) is -a first approximation to either, the -primary^, 
signal or the secondary signal. This first-stage estimacicn 
error signal ej^;^.{t) is input to the first-stage of the 

10 regression filter 80. The first-stage backward prediction 

error b^Ct), multiplied by the first-stage regression, 
coefficient >£:i.x*(t) register 96 value is subtracted from the 
first-stage estimation error signal ei^j,4{t) to produce the 
second-stage estimation error €3^x4 (t:). The second-stage 

15 estimation error signal ej^xa^^) is a second, somewhat better 

approximation to either the primary signal Sx^Ct} or the 
secondary signal nxa(t). 

The same processes are repeated in the least -squares 
lattice predictor 70 and the regression filter 80 for each 

20 stage until a good approximation e^^j,^(t) , to either ' the 

primary signal s^^it) or the secondary signal nj^{t) is 
determined. Each of the signals discussed above, including 
the forward prediction error f„(t) , the backward prediction 
error b„(t) , the estimation error signal e„^x*(t) , is necessary 

25 to calculate the forward reflection coefficient r£^^(t) , the 

backward reflection coefficient rb,„,{t) , and the regression 
coefficient fc^^j^^it) register 90, 92, and 96 values in each 
stage, m. In addition to the forward predictibh' error f^{t), 
the backward prediction error b„(t) , and the estimation error 

30 Smxa(t) signals, a number of intermediate variables, not- 

shown in Figure 8 but based on the values labeled in Figure 
8, are required to calculate the forward reflection 
coefficient r£^„(t), the backward reflection coefficient 
^"b,m(t) , and the regression coefficient J^™,xa(t) register 90, 92, 

35 and 96 values. 
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Intermediate variables include a weighted sum of the 
forward prediction error squares ?>„(t) , a weighted sum of the 
backward prediction error squares S^Ct) , a scalar parameter 
A^{t) , a conversion factor y^{t) , and another scalar parameter 
PBi.x*(t) • The weighted sum of the forward prediction errors 
&aj(t) is defined as: 



a«(t) = i \^'' \fji) (44) 



10 where X without a wavelength identifier^ a or b, is ^ 

constant multiplicative value unrelated to wavelength and is 
typically less than or equal to one, i,e,, X s 1, The- 
weighted sum of the backward prediction errors EmCt) is 
defined as: 

15 

S,(t:) = I X*^'^ |bji) P (45) 



where, again, X without a wavelength identifier, a or b, is 
a constant multiplicative value unrelated to wavelength and 

20 is typically less than or equal to one, i.e., X s 1* These 

weighted sum intermediate error signals can be manipulated, 
such that they are more easily solved for, as described in 
Chapter 9, § 9.3 of the Haykin book referenced above and 
defined hereinafter in equations (59) and (SO) . 

25 DESCRIPTION OF THE JOINT PROCESS ESTIMATOR 

The operation of the joint process estimator 60 is as 
follows. When the joint process estimator 60 is turned on, 
the initial values of intermediate variables and signals 
including the parameter A^.^(t) , the weighted sum of the 

30 forward prediction error signals O^.i(t), the weighted sum of 

the backward prediction error signals Bm.i(t) , the parameter 
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10 



Pin,xa(t)' the zero-Stage estimation error eQ^^it) are 
initialized, some to zero and some to a small positive number 
6: 

A«-i(0) - 0; (46) 

5^i(0) = 6; (47) 

fi„,i(0) - 5; (48) 

P«,x»(0) = 0; (49) 

eoAa^t) = Sxa(t) for t a 0. (50) 

15 After initialization , a simu Itaneoxis sample of the 

measured signal Sj^^lt) or 'S^Cc) and either " rae' secDrfdHry-- 
reference n' (t) or the primary reference s' (t) are input to 
the joint process estimator 60, as shown in Figure 8. The 
forward and backward prediction error signals fo(t) and bo(t) , 

20 and intermediate variables including the weighted sums of the 

forward and backward error signals Oo(t) and So(t), and the 
conversion factor Yo^^) ^.re calculated for the zerjo-st^ge 
according to: 

25 fo(t) = bo(t) r= n' (t) (51a) 

Oo(t) - Eo(t) = X Oo(t-l) + !n'(t)|^ {52a) 

* ♦ 

joit-l) = 1 i53a-> 

30 

if a secondary reference n' (t) is used or according to: 

f,(t) = bo(t) = s' (t) (51b) 

35 Do(t) = /?o(t) = X Oo{t-l) + !s'(t)P (52b) 
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= 1 



(53b) 



if a primary reference s' (t) is used where, again, X without 
a wavelength identifier, a or b, is a constant multiplicative 
5 value unrelated to wavelength- 

Forward reflection coefficient T^^^(t) , backward 
reflection coefficient rb.a{t) , and regression coefficient 
J^a.xjiCt) register 90, 92 and 96 values in .each, f^i-^rj^ 
thereafter are set according to the output of trhe p revious - 

10 stage. The forward reflection coefficient Tf^^ix) , b ackwa rd" 
reflection coefficient Tb^iCt), and regression coefficient 
Ki.xa(t) register 90, 92 and 96 values in the first stage are 
thus set according to the algorithm using values in the zero- 
stage of the joinu process estimator 60. Iti each sLaye, mai", 

15 intermediate values and register values including the 

parameter A x\t} j the forward reflection coe^^^ c*' e*^*" ^ /t-x 
register 90 value; the backward reflection coefficient r^^aCt) 
register 92 value; the forward and backward error signals 
f^it) and b„(t) ; the weighted sum of squared forward 

20 prediction errors 0;^„{t) , as manipulated in § 9.3 of the 

Haykin book; the weighted sum of squared backward prediction 
errors B^^^^it), as manipulated in § 9,3 of the Haykin book; 
the conversion factor 7„(t); the parameter P„,xa(t:); the 
regression coefficient register 96 value; and -the 

25 estimation error e^^ua^^) value are set according to: 

A,.,(t) = X A„.,(t-1) + {b,-i(t-l)fVi(t)/Y„.,(t-l)} (54) 

r,..{t) = -{A„.x(t)/S,.,(t-1)} {-S5v 

30 

r.^Jt) = -{AVi(t)/9f..,(t)} (56) 

fjt) = f..,{t:) + r,^Jt)b,.,(t-l) (57) 

35 b,{t) = b„.,(t-l) + rVm(t)f„.i{t) (58) 
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^Jt) = Vi(t) - {|A».i(t) |VS,.,(t-l)} (59) 

S»(t) = S,.,{t.l) - {|A,,i{t) jV5>„.x(t)} (60) 

S Y„(t-1) - Y..i(t-1) - {|b,.i(t-l) |Vfi«.,(t:-l)} (61) 

P«,xa(t) = Xp„,;^(t-1) + {b,{t)eVx*(t)/7«(t)} (62) 

^,x*(t) = {P«,x*(t)/S«(t)} (^) 

10 

€«.i.x*(t:) = €.^x*(t) - ^•„(t)b„{t) (64) 

where a (*) denotes a complex conjugate* 

These equations ca-uee th-e -error sigrr6rl-s f„^(1r) , 'b-^vtr), 

15 ^tn,xa(t) to be squared or to be multiplied by one another, in 

effect squaring the errors, and r;yg>at|ng np> w i n^P>-rTn^r^^ 
error values, such as A^.^Ct). The error signals and the 
intermediate error values are recursively tied together, as 
shown in the above equations (54) through (64) • They 

20 interact to minimize the error signals in the next stage 

After a good approximation to either the primary signal 
Sxa(t) or the secondary signal n;^a(t) has been determined by 
the joint process estimator 60, a next set of samples, 
including a sample of the measured signal S;^^(t) and a sample. 

25 of either the secondary reference n' (t) or the primary 

reference s' (t) , are input to the joint process estimator 60, 
The re-initialization process does not re-occur, such thaJ^ 
the forward and backward reflection coefficierrt ,^{-t:) arrd'- 
rb^«(t) register 90, 92 values and the regression coefficient 

30 *^m,xa(t) register 96 value reflect the multiplicative values 

required to estimate either the primary signal portion Sy^^it) 
or the secondary signal portion nxa(t) of the sample of Sxavt) 
input previously. Thus, information from previous samples is 
used to estimate either the primary or secondary signal 

35 portion of a present set of samples in each stage. 
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In a more numerically stable and preferred embodiment of 
the above described joint process estimator, a normalized 
joint process estimator is used. This version of the joint 
process estimator normalizes several variables of the above - 
described joint process estimator such that the normalized 
variables fall between -1 and 1. The derivation of the 
normalized joint process estimator is motivated in the Haykin 
text as problem 12 on page 640 by redefining. tbfi= va-7-:-fT^,^h^ 
defined according to the following conditions: 



bit) = 



A,(t) 
V^(t)/3„(t-l) 



This transformation allows the conversion of Equations 
(54) - (64) to the following normalized equations: 

K.. ( t) = V, ( t-i ) [1 - 1 f ( t) 1 ^l^'^i - 1 ( t-i ) I -b,., ( t-i) 1^., ( t) 

h [b,.,(t:-l)-V,(t)f, .,(£)] 

" [i-iv.(t:)!r ii-i^„-.(t:)iT; 

\f..,it:)-l„.,{t)h„.,U-l)] 



[l-|A„-,(t)Mp [l-|b„-.(t-l)M] 

=[l-|V,{c) pli8„.,(t-l) 
T„(t)=Y„.,(t)[l-|b„.,(t) l^j 
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P„(t)=X. 



V„(t)^.(t-i) 



Tl/2 



p„(C-l) +Jb„{t)£,(t) 



Initialization of Normalized Joint Process Estimator 

Let N(t) be defined, as the reference noise input at time 
5 index n and U(t) be defined as combined signal plus noise 

input at time index t the following equations apply (see- 

Haykin, p. 619) : 

1. To initialize the algorithm, at time t=0 set 

^..i{0)=6=10-« 

.0 



2. At each instant til, generate the various zeroth- 
order variables as follows: 

7o(t-l)=l 
/3o(t) =X/3,(t-l) +WMt) 

bo(t)=fo(t)=-^iS= 



3. For regression filtering, initialize the algorithm 
by setting at time index t=0 
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P,(0)=0 

4. At each instant thl, generate the zeroth-order 
variable 

£o(t) =U{t) . 

Accordingly, a normalized joint process estimator can be 
5 used for a more stable system. 

In yet another embodiment, the correlation cancellation 
is performed with a QRD algorithm as shown diagrammatically 
in Figure 8a and as described in detail in Chapter 18 of 
Adaptive Filter Theory by Simon Haykin, published by 
10 Prentice-Hall, -oopyrigiTt 138-6. 

The following equations adapted from the Haykin book 
correspond tx>-the- QPn-TiqTi diagram nf Firpirf> Rn <a}^n adapfpd 
from the Haykin book) . 
Computations 

15 a. Piredictions ; For time t = 1, 2, and 

prediction order m = 1, 2, . . . , M, where M is the final 
prediction order, compute: 
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0,.,(t-l) =X/3„.,(t-2) * |€^.„-x(t-l) P 

e^..(t) = c^.^x(t-l) t) - s;..-i(t-l)X^/'irI^r(-t-D 

irr.,-i(t) = c^...i(t-l)X^'^7r;.„.i(t-l> + s^,„.,(t-l)6^,„.,(t) 
Yr'(t-l) =c,.,.;,(t-l)7i''! (t-l) 



(t) 



^-V^ ( t) 

ei,..(t)=c,,,.i{t)e^,..,{t-i) - s;,,.i(t)x^/=»TT;,,.x{t-i) 

<»-i(t) = c,,„.,{t)X^/=^7r;..-x(t-l)+s,,„.,(t)e^.„.,(t-l) 

b. Filtering : For order m = 0, 1, ... , M-1; and time 
t=l, 2, . . . / compute 

/3„(t) =X^„{C-l)Heb.»(C) P 
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x^/^^y (t-1) 



c. At) = 



p; ^b.- ( t ) x^/'p; ( t-i) +s^., ( t) €. { t) 



Initialization 

a. Auxiliary parameter ini tialization; for order m = 
, M, set 

7r5.,-x(0)=7r^.„.i(0)=0 
p„{0)= 0 

b. Soft copstraint initialization ; For order in = 0, 
. . , M, set • • • 

0J-1)=6 



where 5 is a small positive constant. 

c. Data initialization ; For t = 1, 2, compute 
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€,{t)=d(£:) 
To(t) =1 

where fiit^ is ^Vt^ input and d{t) is the des-ired response 
at time t. 

FLOWCHART OF JOINT PROCESS ESTIMKTOR 
5 In a signal ptrocessor, such as a physiological monitor 

incorporating a reference processor of the present invention 
to determine a reference n' (t) or s' (t) for input to a 
correlation canceler, a joint process estimator 6ti type 
adaptive noise canceler is generally implemented via a 
10 software program havlny an i Lei ati-ve loop. -CSne— iteration -of 
the loop is analogous to a single stage of the joint process 
estimator as shown in Figure 8. Thus^ if a loop is iterated 
m times, it is equivalent to an m stage joint process 
estimator 60 • 

15 A flow chart of a subroutine to estimate the primary 

signal portion s^^it) or the secondary signal portion nxa(t) 
of a measured composite signal, S^aCt) is shown in Figure 9, 
The flow chart illustrates the function of a reference 
processor for determining either the secondary reference 

20 n' (t) or the primary reference s' (t) . The flowchart for the 

joint process estimator is implemented in software. 

A one-time initialization is, performed" when tne 
physiological monitor is powered-on, a^ indicated by sxr 
"INITIALIZE NOISE CANCELER" action block 120. The 

25 initialization sets all registers 90, 92, and 96 and delay 

element variables 110 to the values described above in 
equations (46) through (50) , 

Next, a set of simultaneous samples of the composite 
measured signals S;,^{t) and S^tCt) is input to the subroutine 

30 represented by the flowchart in Figure 9. Then a time update 
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of each of the delay element program variables occurs, as 
indicated in a "TIME UPDATE OF [Z*^] ELEMENTS" action block 
130. The value stored in each of the delay element variables 
110 is set to the value at the input of the delay element 
5 variable 110- Thus, the zero-stage backward prediction error 

bo(t) is stored as the first-stage delay element variable, 
and the first-stage backward prediction error h^it) is stored, 
as the second- stage delay element variable, and so on. 

Then, using the set of measured signal sampi e?=; ^;^(JD ^aniL 

10 Sxb(t) , the reference signal is obtained using the ratiometric 

or the constant saturation methods described above. This i-s" 
indicated by a "CALCULATE REFERENCE [n' (t) or s' (t)] FUR TWO 
MEASURED SIGNAL SAMPLES" action block 140. 

A zero- stage- order update is perfo^r rr'^d r^"^^ mr^^ r^s¥-Ar\ 

15 in a "ZERO-STAGE UPDATE" action block 150. The zero-stage 

backward prediction error bn(t)^ and the zero-stage forward 
prediction error fo(t) are set equal to the value of the 
reference signal n' (t) or s' (t) . Additionally, the weighted 
sum of the forward prediction errors 9f«(t) and the weighted 

20 sum of backward prediction errors fi„(t) are set equal to the 

value defined in equations (47) and (48) . 

Next, a loop counter, m, is initialized as indicared in* 
a "m=0" action block 160. A maximum value of m, defining the 
total number of stages to be used by the subroutine 

25 corresponding to the flowchart in Figure 9, is also defined. 

Typically, the loop is constructed such that it stops 
iterating once a criterion for convergence upon a best 
approximation to either the primary signal ox, the seco^riary 
signal has been met by the joint process estimator 60. 

3 0 Additionally, a maximum number of loop iterations may be 

chosen at which the loop stops iteration. In a preferred 
embodiment of a physiological monitor of the present 
invention, a maximum number of iterations, m=:6 to m=10, is 
advantageously chosen, 

3 5 Within the loop, the forward and backward reflection 

coefficient r^^^it) and r^.^Ct) register 90 and 92 values in 
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the least-squares lattice filter are calculated first, as 
indicated by the "ORDER UPDATE MTH CELL OF LSL-LATTICE" 
action block 170 in Figure 9- This requires calculation of 
intermediate variable and signal values used in determining 
5 register 90, 92, and 96 values in the present stage, the next 

stage, and in the regression filter 80. 

The calculation of regression filter register 96 value 
is performed next, indicated by the "ORDER UPDATE_J£TiL 
STAGE OF REGRESSION FILTER (S ) « action block 180. The two 

10 order update action blocks 170 and 180 are performed in 
sequence m times, until m has reached its predetermined- 
maximum (in the preferred embodiment, m=6 to m=10) or a 
solution has been converged upon, as indicated by a YES path 
from a "DONE" decision block 190, In a computer subroutine, 

15 convergence is determined by checking if the weighted sums of 

the forward and'backward prediction errors tr) and-S-(1r)- are- 
less than a small positive number. An output is calculated 
next, as indicated by a "CALCULATE OUTPUT" action block 200, 
The output is a good approximation to either the primary 

20 signal or secondary signal, as determined by the reference 

processor 26 and joint process estimator 60 subroutine- 
corresponding to the flow chart of Figure 9. This is 
displayed (or used in a calculation in another subroutine) , 
as indicated by a "TO DISPLAY" action block 210. 

25 A new set of samples of the two measured signals 

S;^^(t) and Sxb(t) is input to the processor and joint process 
estimator 60 adaptive noise canceler subroutine corresponding 
to the flowchart of Figure 9 and the process reiterates for 
these samples. Note, however, that the initialization 

3 0 process does not re-occur- New sets of measured signaJ^ 

samples S;^a{t) and Sxb(t) are continuously input to the 
reference processor 26 and joint process estimator adaptive 
noise canceler subroutine. The output forms a chain of 
samples which is representative of a continuous wave. This 

35 waveform is a good approximation to either the primary 

signal waveform S;^^(t) or the secondary wavef orm * n;^, (t) at 
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wavelength Xa . The waveform may also be a good approximation 
to either the primary signal waveform Sxb(t) or the secondary 
waveform n"xb{t) at wavelength Xb. 

A corresponding flowchart for the QRD algorithm of 
Figure 8a is depicted in Figure 9a, with reference numeral 
corresponding in number with an 'a' extension 
CALCULATION OF SATUR ATION FROM 
CORRELATION CANCELER OUTPUT 
Physiological monitors may use the approximation of the- 
primary signals s\Jt) or s"xb(t) or the secondary signals 
n"x,(t) or n"xb<t) to calculate another quantity, such as the 
saturation of one constituent in a volume containing that 
constituent plus one or more other constituents. Generally, 
such calculations require information about either a primary 
or secondary signal at two wavelengths. For example, the 
constant saLuxaLluii metho d-requires a good -approx-xma-ticn-c*. 
the primary signal portions Sx,(t) and s^^{t) of both measured 
signals Syjt) and S^bCt) . The arterial saturation is 
determined from the approximations to both signals, i.e. 
s'\^{t) and s"xb(t) . The constant saturation method also 
requires a good approximation of the secondary signal 
portions nxa(t) or n,b(t) . An estimate of the venous 
saturation may be determined from the approximations to these 
signals i. e. n"x.(t) and n"xb(t) . 

A joint process estimator 60 having two regression 
filters 80a and 80b is shown in Figure 10. A first 
regression filter 80a accepts a measured signal S^^(t) . A 
second regression filter 80b accepts a measur'ed .signal Sxb(c) 
for a use of the constant saturation method to determine th€r 
reference signal n' (t) or s' (t) . The first and secorot 
regression filters 80a and 80b are independent. The backward 
prediction error b„(t) is input to each regression filter 6&a- 
and 80b, the input for the second regression filter 80b 
bypassing the first regression filter 80a. 

The second regression filter 80b comprises registers 98, 
and summing elements 108 arranged similarly to those in the 
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first regression filter 80a. The second regression filter 
80b operates via an additional inteirmediate variable in 
conjunction with those defined by equations (54) through 
(64) , i.e.: 

5 

P«,xb(t:} = Xp^^xbtt-l) + {b„(t)eVxb(t)/7„(t)}; and (65) 

Po,xb(0) = 0. . (66) 

10 The second regression filter 80b has an error signal value 

defined similar to the first regression filter error signal 
values, e„,i^x*(t), i.e.: 



15 



e„,.i,xb(t:) = e„,xb(t) " K\^^{t)hjt) ; and (t?'7} 



^o,xb(t) = Sxb(c"/"i?or-ir a G. (£S) 



The second regression filter has a regression coefficient 
>^m,xb{t) register 98 value defined similarly to the first 
20 regression filter error signal values, i.e.: 

>^«.Ab{t) = {P«,xb(t)/gjt)}; or (69) 

These values are used in conjunction with those intermediate 
25 variable values, signal values, register and register values 
defined in equations (46) through (64) . These signals are 
calculated in an order defined by placing the aniditicnal 
signals immediately adjacent a similar signal for che 
wavelength Xa. 

30 For the constant saturation method, Sxb(t) is input ~t~a 

the second regression filter 80b. The output is then a goorL 
approximation to the primary signal s"xb(^) or secondary"" 
signal s*\^{t) . 

The addition of the second regression filter sab does 

35 not substantially change the computer program subroutine 

represented by the flowchart of Figure 9. Instead of an 
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order update of the m^*^ stage of only one regression filter, 
an order update of the m^^ stage of both regression filters 
80a and 80b is performed. This is characterized by the 
plural designation in the "ORDER UPDATE OF m''*' STAGE OF 
5 REGRESSION FILTER (S) " activity block 180 in Figure 9. Since 

the regression filters 80a and 80b operate independently, 
independent calculations can be performed in the reference 
processor and joint process estimator 60 adaptive -noise- 
canceler subroutine modeled by the flowchart afr F i g ure S . 

10 An alternative diagram for the joint process estimator 

of Figure 10, using the QRD algorithm and having two 
regression filters is shown in Figure 10a. This type of 
joint process estimator would be used for correlation 
cancellation using the QRD algorithm descri'bea m the Hayxxn 

15 book , 

CALCULATION OF -SATURATION 
Once good approximations to the primary signal portions 
s'\^{t) and s«xb(t) secondary signal portion, n"x4{t) 

and n"xb(^)/ have been determined by the joint process 
20 estimator 60, the saturation of Ag in a volume containing. Ag 

and Ag, for example, may be calculated according to various- 
known methods. Mathematically, the approximations to the 
primary signals can be written in terms of Xa and Xb, as: 

25 s\Jt) - es^xaCs^s.^tt) + eg^xaCgXs^sCt) ; and (70) 

s"xb(t:) - €s,xbCsXs,€(t) + fg^xbCfiXs^gCt) , (71) 

Equations (70) and (71) are equivalent to two equationo' 
30 having three unknowns, namely C5(t), C6(t) and . Tixe: 

saturation can be determined by acquiring approximations -to 
the primary or secondary signal portions at two differ ent, 
yet proximate times t^ and tj over which the saturation of Aj 
in the volume containing As and A^ and the saturation of A. 
3 5 in the volume containing A^ and A^ does not change 
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substantially. For example, for the primary signals 
estimated at times t^ and tj: 

s%,(ti) - €s.x.CsXs,s (tj + es,x.C5Xs,e(ti) (72) 

5 

s"xb(ti) - es.xbCsXs.«(ti) + €s.XbC«Xs,«(tJ (73) 

s"x.(t2) - €s^CsXs.j(tj) + e«.x.C(Xs,«(t2) (74) 

Then, difference signals may be determined which relate the 
signals of equations (72) through (75), i.e.: 

15 ASx. = s'\,(tj -s"x,(tj) - €s,x»CsAx + es,x*C6Ax; and (76) 

ASxb = S"xb(ti) -S"xb(t2) - Cs.xbCs^X + e.xbCsAx; (77) 

where Ax = Xs^s(ti) - Xs,e{t2). The average saturation at time 
20 t= (ti+tj) /2 is: 

Saturation(t) = Cs(t)/[Cs(t) + CgCt)]' l"7o; 

€6.X.-€«.Xb{ASx,/ASxb) 

25 = (79) 

Ss.x.-Cs,xa- (€s,xb-es,xb) (ASx»/ASxb) 

It will be understood that the Ax term drops out from the 
saturation calculation because of the division, Thus^ 
3 0 knowledge of the thickness of the primary constituents is not 

required to calculate saturation, 

PULSE OXIMETRY MEASUREMENTS 
A specific example of a physiological monitor utilizing 
a processor of the present invention to' determine a secondary 
35 reference n' (t) for input to a correlation canceler that 

removes erratic motion- induced secondary signal portions is 
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a pulse oximeter. Pulse oximetry may also be performed 
utilizing a processor of the present invention to determine 
a primary signal reference s' (t) which may be used for 
display purposes or for input to a correlation canceler to 
5 derive information about patient movement and venous blood 

oxygen saturation. 

A pulse oximeter typically causes energy to propagate 
through a medium where blood flows close to the surface for 
example, an ear lobe, or a digit such as a f^ngpr, r fo>-^h*>Rrl 
10 or a fetus' scalp. An attenuated signal is measured aft-er 
propagation through or reflected from tTie medium. The pulse 
oximeter estimates the saturation of oxygenated blood. 

Freshly oxygenated blood is pumped at high pressure fzom 
• the heart into the arteries for use by the bcd^'-, Ihs_-zcl™_ 
15 of blood in the arteries varies with the heartbeat, giving- 

xise to a variation in absorption of enercry. at the rate of 
the heartbeat, or the pulse. 

Oxygen depleted, or deoxygenated, blood is returned to 
the heart by the veins along with unused oxygenated blood. 
20 The volume of blood in the veins varies with the rat.e of 

breathing, which is typically much slower than the heartbeat - 
Thus, when there is no motion induced variation in: Llic 
thickness of the veins, venous blood causes a low frequency 
variation in absorption of energy. When there is motion 
.25 induced variation in the thickness of the veins, the low 

frequency variation in absorption is coupled with the erratic 
variation in absorption due to motion artifact. 

In absorption measurements using the traxasmi ssirin ri£ 
energy through a medium, two light emitting diodes (LED's) 
3 0 are positioned on one side of a portion oL the body where 

blood flows close to the surface, such as a finger, and a 
photodetector is positioned on the opposite side q£ the 
finger. Typically, in pulse oximetry measurements, one LED 
emits a visible wavelength, preferably red, and the other LED 
35. emits an infrared wavelength. However, one skilled in the 
art will realize that other wavelength combinations could be 
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used. The finger comprises skin, tissue, muscle, both 

arterial blood and venous blood, fat, etc., each of which 
absorbs light energy differently due to different absorption 
coefficients, different concentrations, different 
5 thicknesses, and changing optical pathlengths. When the 

patient is not moving, absorption is substantially constant 
except for the flow of blood. The constant attenuation can 
be determined and subtracted from the signal via tradi±±Qnal. 
filtering techniques. When the patient moves, this ca useb 

10 perturbation such as changing optical pathlength due to' 

movement of background fluids (e,g., venous blood having a 
different saturation than the arterial blood) . Therefore, 
the measured signal becomes erratic. Erratic motion induced 
noise typically cannot be predetermined and/or sub uiacLed ' 

15 from the measured signal via traditional filtering 

techniques* Thtrs-, determining the- cxj'^gen — sarturaticn -of 
arterial blood and venous blood becomes more difficult. 

A schematic of a physiological monitor for pulse 
oximetry is shown in Figures 11-13. Figure 11 depicts a 

20 general hardware block diagram of a pulse oximeter 299-. . A 

sensor 300 has two light emitters 301 and 302 such as LSD's. 
One LED 301 emitting light of red wavelengths and another LED 
302 emitting light of infrared wavelengths are placed 
adjacent a finger 310. A photodetector 320, which produces 

25 an electrical signal corresponding to the attenuated visible 

and infrared light energy signals is located opposite the 
LED's 301 and 302. The photodetector 320 is connected to 
front end analog signal conditioning circuity *330 . 

The front end analog signal conditioning circuitry 33G 

30 has outputs coupled to analog to digital conversion circuit 

332. The analog to digital conversion circuitry 332 has 
outputs coupled to a digital signal processing system 334. 
The digital signal processing system 334 provides the desired 
parameters as outputs for a display 336. Outputs for'display 

35 are, for example, blood oxygen saturation, heart rate, and a 

clean plethysmographic waveform. 
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The signal processing system also provides an emitter 
current control output 337 to a digital-to-analog converter 
circuit 338 which provides control information for light 
emitter drivers 340. The light emitter drivers 340 couple to 
5 the light emitters 301, 302, The digital signal processing 

system 334 also provides a gain control output 342 for the 
front end analog signal conditioning circuitry 330. 

Figure 11a illustrates a preferred embodiment for the 
combination of the emitter drivers 340 and the digital to 

10 analog conversion circuit 338* As depicted in Figure 11a, 

the driver comprises First and second input latches 321, 322, 
a synchronizing latch 323, a voltage reference 324, a digital 
to analog conversion circuit 325, first and second switch 
banks 326, 327, firet and second voltage- to current 

15 converters 328, 329 and the LED emitters 3 01, 302 

corresponding to the LED emitters 301, 302 of Figure _1JL.. 

The preferred driver depicted in Figure 11a is 
advantageous in that the present inventors recognized that 
much of the noise in the oximeter 299 of Figure 11 is caused 

20 by the LED emitters 301, 302. Therefore, the emitter driver 

circuit of Figure 11a is designed to minimize the noise from 
the emitters 301, 302, The first and second input Tatches 
321, 324 are connected directly to the DSP bus • Therefore, 
these latches significantly minimizes the bandwidth 

25 (resulting in noise) present on the DSP bus which passes 

through to the driver circuitry of Figure 11a, The output of 
the first and second input latches only changes when these 
latched detect their address on the DSP bus. The , first input 
latch receives the setting for the digital La analog. 

30 converter circuit 325. The second input latch receives 

switching control data for the switch banks 326, 327. The 
synchronizing latch accepts the synchronizing pulses which 
maintain synchronization between the activation of emitters 
301, 302 and the analog to digital conversion circuit 332. 

3 5 I'he voltage reference is also chosen as a low noise DC 

voltage reference for the digital to analog conversion 
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circuit 325, In addition, in the present embodiment, the 
voltage reference has an lowpass output filter with a very 
low corner frequency {e,g., IHz in the present embodiment). 
The digital to analog converter 325 also has a lowpass filter 
5 at its output with a very low coimer frequency (e.g.^ 1 Hz) , 
The digital to analog converter provides signals for each of 
the_ emitters 301, 302. 

In the present embodiment, the output of the voltage to 
current converters 328, 329 are switched such- that- with ■ f.hfi .. 
10 emitters 301, 302 connected in back-to-back c o nf i gur a t i e n r - 
only one emitter is active an any given time. Tn addition, 
the voltage to current converter for the inactive emitter is 
switched off at its input as well, such that it is completreiy 

deactivated- TYl 1 ° T^^WC^^ noise f v' OTn t-Vxa cwt »-r^Vii nr r anH 

15 voltage to current conversion circuitry. In tire pres eii t 

embodiment, low noise voltage to current converters are 
selected (e-g-. Op 27 Op Amps), and the feedback loop is 
configured to have a low pass filter to reduce noise. In the 
present embodiment, the low pass filtering function of the 

20 voltage to current converter 328, 329 has a corner frequency 

of just above 625 Hz, which is the switching speed for the 
emitters, as further discussed below. Accordingly", izhe 
preferred driver circuit of Figure 11a, minimizes the noise 
of the emitters 301, 302, 

25 In general, the red and infrared light emitters 3 01, 3 02 

each emits energy which is absorbed by the finger 310 and 
received by the photodetector 320. The photodetector 320 
produces an electrical signal which corresponds. ±ja th^ 
intensity of the light energy striking the photodetector 3^0. 

30 The front end analog signal conditioning circuitry 330 

receives the intensity signals and filters and conditions 
these signals as further described below fox iuxther 
processing. The resultant signals are provided to the 
analog- to-digital conversion circuitry 332 which converts the— 

35 analog signals to digital signals for further processing by 

the digital signal processing system 334, The digital signal 
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processing system 334 utilizes the two signals in order to 
provide a what will be called herein a "saturation 
transform." It should be understood, that for parameters 
other than blood saturation monitoring, the saturation 
5 transform could be better termed as a concentration 

transform, in-vivo transform, or the like, depending on the 
desired parameter. The term saturation transform is used to 
describe an operation which converts the sample data from 
time domain to saturation domain values as will be ap^iaziant- 
10 from the discussion below. In the present embodiment, -the- 

output of the digital signal processing system 334 provides 
clean plethysmographic waveforms of the detected signals and 
provides values for oxygen saturation and pulse rate tro uhe- 
display 336. 

15 It should be understood that in different embodiments -cTf" 

the present invention, one or more of the outputs may be 
provided. The digital signal processing system 334 also 
provides control for driving the light emitters 301, 3 02 with 
an emitter current control signal on the emitter current 

20 control output 337. This value is a digital value which. is 

converted by the digital-to-analog conversion circuit 338 
which provides a control signal to the emictier curieuu 
drivers 340. The emitter current drivers 340 provide the 
appropriate current drive for the red emitter 301 and the 

25 infrared emitter 302. Further detail of the operation of the 

physiological monitor for pulse oximetry is explained below. 

In the present embodiment, the light emitters are driven 
via the emitter current ' driver 340 to provide light 
transmission with digital modulation at 625 Hz:. In — the 

3 0 present embodiment, the light emitters 3 01, 3 02 are driven at 

a power level which provides an acceptable intensity for 
detection by the detector and for conditioning ixy the- front, 
end analog signal conditioning circuitry 330. Once this 
energy level is determined for a given patient by the digital 

35 signal processing system 334, the current level for the red 

and infrared emitters is maintained constant. It should be 
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understood, however, that the current could be adjusted for 
changes in the ambient room light and other changes which 
would effect the voltage input to the front end analog signal 
conditioning circuitry 330- In the present invention, the 
5 red and infrared light emitters are modulated as follows: for 
one complete 625 Hz red cycle/ the red emitter 301 is 
activated for the first quarter cycle, and off for the 
remaining three-quarters cycle; for one comglete £25 Hz- 
infrared cycle ^ the infrared light emitter 50z -±s dCLlVdL i bfd - 

10 for one quarter cycle and is off for the remaining three- 
quarters cycle; In order to only receive one signal at a 
time, the emitters are cycled on and off alternatively, in 
sequence, with each only active for a quarter cycle per 625 
Hz cycle and a quarter cycle separating the a'crive rimes. 

15 The light signal is attenuated (amplitude modulated) by 

medium) • The attenuated (amplitude modulated) signal is 
detected by the photodetector 320 at the 625 Hz carrier 
frequency for the red and infrared light. Because only a 

20 single photodetector is used, the photodetector 320 receives 

both the red and infrared signals to form a composite time- 
division signal. 

The composite time division signal is provided to the 
front analog signal conditioning circuitry 330. Additional 

25 detail regarding the front end analog signal conditioning 

circuitry 330 and the analog to digital converter circuit 332 
is illustrated in Figure 12. As depicted in Figure 12, the 
front end circuity 302 has a preamplifier T^, - a high pass 
filter 344, an amplifier 346, a programmable gain amplifier 

30 348, and a low pass filter 350, The preamplifier 342 is* 

transimpedance amplifier that converts the composite current 
signal from the photodetector 320 to a corresponding VT^ltrage- 
signal, and amplifies the signal. In the present embodiment, 
the preamplifier has a predetermined gain to boost the signaT 

3 5 amplitude for ease of processing. In the present embodiment, 

the source voltages for the preamplifier 342 are -15 VDC and 
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+ 15 VDC. As will be understood, the attenuated signal 
contains a component representing ambient light as well as 
the component representing the infrared or the red light as 
the case may be in time* If there is light in the vicinity 
5 of the sensor 300 other than the red and infrared light, this 

ambient light is detected by the photodetector 320. 
Accordingly, the gain of the preamplifier is selected in 
order to prevent the ambient light in the. fij gnal from 
saturating the preamplifier under normal and reasonable- 

10 operating conditions. 

In the present embodiment, the preamplifier 342 
comprises an Analog Devices AD743JR OpAmp. This 
transimpedance amplifier is particularly advantageous din th a t 
it exhibits several desired featrxnres" £01 1rite syotrcm' 

15 described/ such as: low equivalent input voltage noise, low 

input current -noi-se, low- input -bia.s ^current , —hisCfh 

gain bandwidth product, low total harmonic distortion, high 
common mode rejection, high open loop gain, and a high power 
supply rejection ratio • 

20 The output of the preamplifier 342 couples as an input 

to the high pass filter 344, The output of the preamplifier 
also provides a first input 346 to the analog to digital 
conversion circuit 332, In the present embodiment, the high 
pass filter is a single-pole filter with a comer frequency 

25 of about 1/2 - 1 Hz* However, the corner frequency is 

readily raised to about 90 Hz in one embodiment ♦ As will be- 
understood, the 625 Hz carrier frequency of the red and 
infrared signals is well above a 90 Hz corner f r^queacy . The" 
high-pass filter 344 has an output coupled as an input tc an 

30 amplifier 346. In the present embodiment, the amplifier 3-4-6- 

comprises a unity gain amplifier. However, the gain of the 
amplifier 346 is adjustable by the variation of a singl-e- 
resistor. The gain of the amplifier 346 would be increased 
if the gain of the preamplifier 342 is decreased to 

35 compensate for the effects of ambient light. 
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The output of the amplifier 346 provides an input to a 
programmable gain amplifier 348* The programmable gain 
amplifier 348 also accepts a programming input from the 
digital signal processing system 334 on a gain control signal 
5 line 343 • The gain of the programmable gain amplifier 348 is 

digitally programmable. The gain is adjusted dynamically at 
initialization or sensor placement for changes in the medium 
under test from patient to patient. For example, the signal 
from different fingers differs somewhat. Therefore, -a 

10 dynamically adjustable amplifier is provided by the 
programmable gain amplifier 348 in order to obtain a signal 
suitable for processing. 

The programmable gain amplifier is also advantageous in 
an alternative embodimenc in which the tJuuLLex drivfe xnxrrentr 

15 is held constant. In the present embodiment, the emitter 

drive current is-adj-U^feed for ^.ar'V> ,rva^^ ^nh _nn._oyri^r' ^r> r^^.^a^n 

the proper dynamic range at the input of the analog to 
digital conversion circuit 332. However, changing the 
emitter drive current can alter the emitter wavelength, which 

20 in turn affects the end result in oximetry calculations. 

Accordingly, it would be advantageous to fix the emi.tter 
drive current for all patients. In an alternative embodiment 
of the present invention, the programmable gain amplifier can 
be adjusted by the DSP in order to obtain a signal at the 

25 input to the analog to digital conversion circuit which is 

properly within the dynamic range (+ 3v to -3v in the present 
embodiment) of the analog to digital conversion circuit 332, 
In this manner, the emitter drive current could be fixed for 
all patients, eliminating the wavelength shift due to emitter 

30 current drive changes. 

The output of the programmable gain amplifier 348 
couples as an input to a low-pass filter 350. 
Advantageously, the low pass filter 350 is a single-pole 
filter with a corner frequency of approximately 10 Khz in the 

35 present embodiment. This low pass filter provides anti- 

aliasing in the present embodiment. 
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The output of the low-pass filter 350 provides a second 
input 352 to the analog- to-digital' conversion circuit 332. 
Figure 12 also depicts additional defect of the analog-to- 
digital conversion circuit. In the present embodiment, the 
5 analog-to-digital conversion circuit 332 comprises a first 
analog- to-digital converter 354 and a second analog- to- 
digital converter 356. Advantageously, the first analog- to- 
digital converter 354 accepts input from the first input 346 
to the analog-to-digital conversion circuit 332, and th-e- 

10 second analog to digital converter 356 accepts input on the 

second input 352 to the analog-to-digital conversion 
circuitry 332, 

In one advantageous embodiment, the first analog- to- 
digital convex ter 354 i'S a diagnostiic analog— to — d"i*gs;t-al 

15 converter. The diagnostic task (performed by the digital 

signal pr-or-^gc;T ng syfitfim)- i s -to read the ^output rtf ^hf> 

detector as amplified by the preamplifier 342 in order to 
determine if the signal is saturating the input to the high- 
pass filter 344 • In the present embodiment, if the input to 

20 the high pass filter 344 becomes saturated, the front, end 

analog signal conditioning circuits 33 0 provides a *0' 
output. Alternatively, the first analog-to-digital converter 
3 54 remains unused. 

The second analog- to-digital converter 352 accepts the 

25 conditioned composite analog signal from the front end signal 
conditioning circuitry 330 and converts the signal to digital 
form. In the present embodiment, the second analog to 
digital converter 356 comprises a single-channel,, delta-sigma 
converter. , In the present embodiment, a Crystal 

30 Semiconductor CS5317-KS delta-sigma analog to digital 

converter is used. Such a converter is advantageous in that 
it is low cost, and exhibits low noise characteristics. More 
specifically, a delta-sigma converter consists of two major 
portions, a noise modulator and a decimation filter. The 

35 selected converter uses a second order analog delta-sigma 

modulator to provide noise shaping. Noise shaping refers to 
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changing the noise spectrum from a flat response to a 
response where noise at the lower frequencies has been 
reduced by increasing noise at higher frequencies. The 
decimation filter then cuts out the reshaped, higher 
5 frequency noise to provide 16 -bit performance at a lower 

frequency* The present converter samples the data 128 times 
for every 16 bit data word that it produces. In this manner, 
the converter provides excellent noise rejection, dynamic 
range and low harmonic distortion, that help in critic:>l 
10 measurement situations like low perfusion and electrocaiitery* 

In addition, by using a single-channel converter, there 
is no need to tune two or more channels to each other. The 
delta-sigma cornrerrt^rr i'S also ad:vaTrt:agecrtrs in -t***o.t: it 

15 exhibits noise shaping, for improved noise control. An 

exemplary analog txx. digifril rnnverter Is — ^ — Cry-sJial_. 

Semiconductor CS5317. In the present embodiment, the second 
analog to digital converter 356 samples the signal at a 20 
Khz sample rate. The output of the second analog to digital 

20 converter 356 provides data samples at 20 Khz to the digital 

signal processing system 334 {Figure 11) . 

The digital signal processing system 334 is illustrated 
in additional detail in Figure 13. In the present 
embodiment, the digital signal processing system comprises a 

25 microcontroller 360, a digital signal processor 3 62, a 

program memory 364, a sample buffer 366, a data memory 3 68, 
a read only memory 370 and communication registers 372. In 
the present embodiment, the digital signal prorce^ox 352 
an Analog Devices AD 21020, " In the present embcdimexii^, the- 

30 microcontroller 360 comprises a Motorola 68HC05, with built 

in program memory. In the present embodiment, the sample 
buffer 366 is a buffer which accepts the 20 Khz -sa^p.pl-e dat-a- 
from the analog to digital conversion circuit 332 for storage 
in the data memory 3 68. In the present embodiment, rhe datra" 

35 memory 368 comprises 32 KWords (words being 40 bits in the 

present embodiment) of static random access memory. 
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The microcontroller 3 60 is connected to the DSP 3 62 via 
a conventional JTAG Tap line. The microcontroller 360 
transmits the boot loader for the DSP 362 to the program 
memory 364 via the Tap line, and then allows the DSP 362 to 
5 boot from the program memory 364, The boot loader in program 
memory 364 then causes the transfer of the operating 
instructions for the DSP 362 from the read only memory 3:7.0_tc^ 
the program memory 364. Advantageously, the program :^.ory 
364 is a very high speed memory for the DSP 3G2. 
10 The microcontroller 360 provides the emitter current 

control and gain control signals via the communications 
register 372. 

Figures 14-20 depict functional block diagrams of the 
operations ot tne pulse oximeter 299 carried out by tne 

15 digital signal processing system 334* The signal process inq 
functions de s cr ibe d-be-l-ow are cairri-ed-out-by-the DS? 3o2-irr" 
the present embodiment with the microcontroller 360 providing 
system management. In the present embodiment, the operation 
is software/firmware controlled. Figure 14 depicts a 

20 generalized functional block diagram for the operations 

performed on the 20 Khz sample data entering the digital 
signal processing system 334. As illustrated in Figure 14, 
a demodulation, as represented in a demodulation module 400, 
is first performed. Decimation, as represented in a 

25 decimation module 402 is then performed on the resulting 

data. Certain statistics are calculated, as represented in. 
a statistics module 404 and a saturation transform i-s 
performed, as represented in a saturation transtorm module 
406, on the data resulting from the dec i mat ion oy e r a L Ion* 

3 0 The data subjected to the statistics operations and the dai::a- 

subjected to the saturation transform operations are 
forwarded to saturation operations, as represented by a 
saturation calculation module 408 and pulse rate operations, 
as represented in a pulse rate calculation module 410. 

3 5 In general, the demodulation operation separates the red 

and infrared signals from the composite signal and removes 
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the 625 Hz carrier frequency, leaving raw data points. The 
raw data points are provided at 625 Hz intervals to the 
decimation operation which reduces the samples by an order of 
10 to samples at 62,5 Hz. The decimation operation also 
provides some filtering on the samples. The resulting data 
is subjected to statistics and to the saturation transform 
operations in order to calculate a saturation value which is 
very tolerant to motion artifacts and othar. noi sp in the^ 
signal. The saturation value is ascertrctined ±n -tfe 



10 saturation calculation module 408, and a pulse rate and -ar 

clean plethysmographic waveform is obtained through the pulse 
rate module 410. Additional detail regarding the various 
operations is provided in connection with Figures 15-21, 

Figure XB" illuscraxes rhe operation or the de i uouulaLluu 

IS module 400, The modulated signal format is depicted in 

Figure* 15. On-e—full— S25 -Hz cycle of the -composite -signal —i3- 
depicted in Figure 15 with the first quarter cycle being the 
active red light plus ambient light signal, the second 
quarter cycle being an ambient light, signal, the third 

20 quarter cycle being the active infrared plus ambient light 

signal, and the fourth quarter cycle being an ambient light 
signal. As depicted in Figure 15, with a 20 KHz sampling 
frequency, the single full cycle at 625 Hz described above 
comprises 32 samples of 20 KHz data, eight samples relating- 

25 to red plus ambient light, eight samples relating to ambient 

light, eight samples relating to infrared plus ambient lights 
and finally eight samples related to ambient light* 

Because the signal processing system 33:^, comnrol^s" trhe' 
activation of the light emitters 300, 302, the-entir^ sys-t-e^. 

30 is synchronous. The data is synchronously divided (and 

thereby demodulated) into four 8 -sample packets, with a time 
division demultiplexing operation as represented in a 
demultiplexing module 421. One eight-sample packet 422 
represents the red plus ambient light signal; a second eight - 

35 sample packet 424 represents an ambient light signal; a third 

eight -sample packet 426 represents the attenuated infrared 
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light plus ambient light signal; and a fourth eight-sample 
packet 428 represents the ambient light signal. A select 
signal synchronously controls the demultiplexing operation so 
as to divide the time-division multiplexed composite signal 
5 at the input of the demultiplexer 421 into its four subparts. 

A sum of the last four samples from each packet is then 
calculated, as represented in the summing operations 430, 
432, 434, 43 6 of Figxire 15. In the present emhodimpnf . . ±h&. 
last four samples are used because a low -psss filte^r xrr LLe- 
10 analog to digital converter 356 of the present embudiuitiiiL lias 
a settling time. Thus, collecting the, last four samples from 
each 8 -sample packet allows the previous signal to clear. 
This summing operation provides an integration operation 
which enhances noise immunity. The sutir ul Lue- LebyeCLrive- 
15 ambient light samples is then subtracted from the sum of the 

red and inf rered-saciples, -as - represented -in -the -subtrac-t-ion. 
modules 438, 440. The siibtraction operation provides some 
attenuation of the ambient light signal present in the data. 
In the present embodiment, it has been found that 
20 approximately 20 dB attenuation of the ambient light -is 
provided by the operations of the subtraction modules 438, 
440. The resultant red and infrared sum values are divided 
by four, as represented in the divide by four modules 442, 
444. Each resultant value provides one sample each of the 
25 red and infrared signals at 625 Hz. 

It should be understood that the 625 Hz carrier 
frequency has been removed by the demodulation operation 40n 
The 625 Hz sample data at the output ox -trip, demodulatrixnr 
operation 400 is sample data without the ca-rrie-r frequency. 
3 0 In order to satisfy Nyquist sampling requirements, less than 

20 Hz is needed (understanding that the human pulse is about 
25 to 250 beats per minute, ' or about .4 Hz - 4 -Hz) . 
Accordingly, the 625 Hz resolution is reduced to 62.5 Hz in 
the decimation operation. 
3 5 Figure 16 illustrates the operations of the decimation 

module 402. The red and infrared sample data is provided at 
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625 Hz to respective red and infrared buffer/filters 450, 
452. In the present embodiment, the red and infrared 
buffer/filters are 519 samples deep. Advantageously, the 
buffer filters 450, 452 function as continuous first-- in, 
5 first *out buffers* The 519 samples are subjected to low-pass 

filtering. Preferably, the low-pass filtering has a cutoff 
frequency of approximately 7.5 Hz with attenuation of 
approximately -110 dB* The buffer/filters 450, 452 form h 
Finite Impulse Response (FIR) filter with coefficients fui 

10 519 taps. In order to reduce the sample frequency by ten, 

the low-pass filter calculation is performed every ten 
samples, as represented in respective red and infrared 
decimation by 10 modules 454, 456. In other words, with the 
transfer of each new ten samples into the bufTer/'fi"Iters 4"5'0, 

15 452, a new low pass filter calculation is performed by 

multiplying tiie i ' uipulse respons e- (troe-f f i-ci-ent-s) by-the— 
filter taps. Each filter calculation provides one output 
sample for respective red and infrared output buffers 458, 
460. In the present embodiment, the red and infrared output 

20 buffers 458, 460 are also continuous FIFO buffers that hold 

570 samples of data. The 570 samples provide respectiv-e- 
infrared and red samples or packets (also denoted "snapshot" 
herein) of samples. As depicted in Figure 14, the output 
buffers provide sample data for the statistics operation 

25 module 404, saturation transform module 406, and the pulse 

rate module 410. 

Figure 17 illustrates additional functional operation 
details of the statistics module 404. In , summary, che 
statistics module 404 provides first order oximetry" 

30 calculations and RMS signal values for the red and infraredr 

channels. The statistics module also provides a cross- 
correlation output which indicates a cross-correlaticmr 
between the -red and infrared signals. 

As represented in Figure 17, the statistics operation 

35 accepts two packets of samples (e.g., 570 samples at 62.5 Hz 

in the present embodiment) representing the attenuated 
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infrared and red signals, with the carrier frequency removed. 
The respective packets for infrared and red signals are 
normalized with a log function, as represented in the Log 
modules 480, 482 • The normalization is followed by removal 
5 of the DC portion of the signals, as represented in the DC 

Removal modules 484, 486. In the present embodiment, the DC 
removal involves ascertaining the DC value of the first one 
of the samples (or the mean of the first several or the mean, 
of an entire snapshot) from each of the respective red and 

10 infrared snapshots, and removing this DC value from all 

samples in the respective packets. 

Once the DC signal is removed, the signals are subjected 
to bandpass filtering, as represented in red and infrared. 
Bandpass Fil Lex ruodul ca^^'oS, 430. In t^re -preKenf cuiLudimciiL, 

15 with 570 samples in each packet, the bandpass filters are 

configured JwitJi^£^-l— tap s t o - provi d e -a^ FXR--iiJ.t^. -Viitii— a- 
linear phase response and little or no distortion. In the 
present embodiment, the bandpass filter has a pass band from 
34 beats/minute to 250 beats/minute. The 301 taps slide over 

20 the 570 samples in order to obtain 270 filtered samples 

representing the filtered red signal and 270 filtered samples 
representing the filtered infrared signal. In an ideal case, 
the bandpass filters 488, 490 remove the DC in the signal. 
However, the DC removal operations 484, 486 assist in DC 

25 removal in the present embodiment. 

After filtering, the last 120 samples from each packet 
(of now 270 samples in the present embodiment) are selected 
for further processing as represented in Sele^ct: Last: 120 
Samples modules 492, 494. The last 120 samples are selected 

30 because, in the present embodiment, the first 150 samples 

fall within the settling time for the Saturation Transfer 
module 406, which processes the same data packets, as furthe~r 
discussed below. 

Conventional saturation equation calculations are 

35 performed on the red and infrared 120-sample packets. In the 

present embodiment, the conventional saturation calculations 
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are performed in two different ways. For one calculation, 
the 120-sample packets are processed to obtain their overall 
RMS value, as represented in the first red and infrared RMS 
modules 496, 498. The resultant RMS values for red and 
infrared signals provide input values to a first 
RED_RMS/IR_RMS ratio operation 500, which provides the RMS 
red value to RMS infrared value ratio as an input to a 
saturation equation module 502. As well understood in the 
art, the ratio of the intensity of red to infrared ^tt^.ua^d 
light as detected for known red and infrared waveiengths- 
(typically = 650 nm and Xjr = 910 nm) relates to the 

oxygen saturation of the patient. Accordingly, the 
saturation equation module 502 represents a conventional 
look-up table or Lhc Ixke which, for- piredet^TTirrricd— Xairxos7 
provides known saturation values at its output 504. The red 
and infrared RMS- v^lnp" r^r e^ a1 q n prnviripfi rs nntpnt-s of ... t^^— 
statistics operations module 404. 

In addition to the conventional saturation operation 
502, the 120-sample packets are subjected to a cross- 
correlation operation as represented in a first cross - 
correlation module 506. The first cross -correlation module. 
506 determines if good correlation exists between the 
infrared and red signals. This cross correlation is 
advantageous for detecting defective or otherwise 
malfunctioning detectors. The cross correlation is also 
advantageous in detecting when the signal model (i.e., the 
model of Equations (1) - (3)) is satisfied. Xf r nrrel a tio n 
becomes too low between the two channels, the sxgrrarl -rnodel xs" 
not met. In order to determine this, the normalized cross, 
correlation can be computed by the cross -correlation, module 
506 for each snapshot of data. One such correlation function 
is as follows: 



If the cross correlation is too low, the oximeter 299 
provides a warning (e.g., audible, visual, etc) to the 
operator. In the present embodiment, if a selected snapshot 
yields a normalized correlation of less than 0.75, the 
5 snapshot does not qualify. Signals which satisfy the signal 
model will have a correlation greater than the threshold. 

The red and infrared 120 -sample packets are also 
subjected to a second saturation operation and cross 
correlation in the same manner as described abovB., f^yrpp!" ^h^ 
10 120 samples are divided into 5 equal bins of samples {i.e., 

5 bins of 24 samples each) . The RMS, ratio, saturation, and" 
cross correlation operations are performed on a bin-by-bin 
basis. These operations are represented in the Divide xnco 

TT" •* • 1^1 D J *«3 • ■» T i— C T A CIO ^ e-fix/^ ^ v% ^ ^"i^^ 4 ^ ^ ^^-^ *^ 

15 RMS modules 514, 516, the second RED-RMS/IR-RMS ratio modure 

518, the second saturation equation module 520 and the second 
cross correlation module 522 in Figure 17. 

Figure 18 illustrates additional detail regarding the 
saturation transform module 406 depicted in Figure 14* As 

20 illustrated in Figure 18, the saturation transform module 406 

comprises a reference processor 530, a correlation canceler 
531, a master power curve module 554, and a bin power curve 
module 533. The saturation transform module 406 can be 
correlated to Figure 7a which has a reference processor 2,6 

25 and a correlation canceler 27 and an integrator 29 to provide 
a power curve for separate signal coefficients as depicted in 
Figure 7c. The saturation transform module 406 obtains a 
saturation spectrum from the snapshots of det-a , In other 
words, the saturation transform 4 06 provides i nf n•rm^ ^ -i o-r 

30 the saturation values present in the snapshots. 

As depicted in Figure 18, the reference processor 530 
for the saturation transform module 406 has a saturation- 
equation module 532, a reference generator module 534, a DC 
removal module 536 and a bandpass filter module 538. T he red " 

35 and infrared 570-sample packets from the decimation operation 

are provided to the reference processor 530. In addition, a 
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plurality of possible saturation values (the "saturation axis 
scan") are provided as input to the saturation reference 
processor 530, In the present embodiment, 117 saturation 
values are provided as the saturation axis scan. In a 
5 preferred embodiment, the 117 saturation values range 

uniformly from a blood oxygen saturation of 34*8 to 105«0. 
Accordingly, in the present embodiment, the 117 saturation 
values provide aa axis scan for the reference prnnessor 530. 
which generates a reference signal for use by t Axe* corr^'latzfotv 

10 canceler 531. In other words, the reference processoir ±s^ 

provided with each of the saturation values, and a resultant 
reference signal is generated corresponding to the saturation 
value- The correlation canceler is formed by a joint process 
estimator 5313 and a low pass filter 552 ±n- trhe pxefserrt*-*" 

15 embodiment. 

It should ^be und^^s-t-ocd that the -scan— values could— be. 

chosen to provide higher or lower resolution than 117 scan 
values. The scan values could also be non-uniformly spaced. 
As illustrated in Figure 18, the saturation equation 

20 module 532 accepts the saturation axis scan values as .an 

input and provides a ratio "r„" as an output. In comparison 
to the general discussion of Figure 7a-7c, this ratio "r^" 
corresponds to the plurality of scan value discussed above in 
general. The saturation equation simply provides a known 

25 ratio "r" (red/infrared) corresponding to the saturation 
value received as an input. 

The ratio "r^" is provided as an input to the reference 
generator 534, as are the red and infraretx sample packecy ; 
The reference generator 534 multiplies either the- red or 

3 0 infrared samples by the ratio "r^^" and subtracts the value 

from the infrared or red samples, respectively. For 
instance, in the present embodiment, the reference generator 
534 multiplies the red samples by the ratio "r^" and 
subtracts- this value from the infrared samples. The 

35 resulting values become the output of the reference generator 

534, This operation is completed for each of the saturation 
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scan values (e.g., 117 possible values in the present 
embodiment) , Accordingly, the resultant data can be 
described as 117 reference signal vectors of 570 data points 
each, hereinafter referred to as the reference signal 
vectors. This data can be stored in an array or the like. 

In other words, assuming that the rfed and infrared 
sample packets represent the red S^edtt) and infrared Sia(t) 
measured signals which have primary s(t) and secondary jaXD. 
signal portions, the output of the reference generoctor- 
becomes the secondary reference signal n' (t) , which compl±HS- 
with the signal model defined above, as follows: 

n'it) =S,,{t) -r„S,^(t) 

In the present embodiment, the refereiice-srigilal vecLUis 
and the infrared signal are provided as input to the DC 
removal module 536-of the reference processor 53Q. The-DC 
removal module 536, like the DC removal modules 484, 486 in 
the statistics module 404, ascertains the DC value of the 
first of the samples for the respective inputs (or mean of 
the first several or all samples in a packet) and subtracts 
the respective DC baseline from the sample values. The 
resulting sample values are subjected to a bandpass filter 
538. 

The bandpass filter 538 of the reference processor 530 
performs the same type of filtering as the bandpass filters 
488, 490 of the statistics module 404. Accordingly, each set 
of 570 samples subjected to bandpass filtering results in 22IL 
remaining samples. The resulting data at a tirsi; output 542 
of the bandpass filter 538 is one vector of 270 samples 
(representing the filtered infrared signal in the present 
embodiment) . The resulting data at a second output 540 of 
the bandpass filter 538, therefore, is 117 reference signal 
vectors of 270 data points each, corresponding to each of the 
saturation axis scan values provided to the saturation 
reference processor 53 0. 
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It should be understood that the red and infrared sample 
packets may be switched in their use in the reference 
processor 530. In addition, it should be understood that the 
DC removal module 536 and the bandpass filter module 538 can 
5 be executed prior to input of the data to the reference 
processor 530 because the calculations performed in the 
reference processor are linear. This results in a 
significant processing economy. 

The outputs of the reference processor 530 provide firsir 

10 and second inputs to a joint process estimator 550 of the 
type described above with reference to Figure 8. The first 
input to the joint process estimator 550 is the 270-sample 
packet representing the infrared signal in the present 
embodiment. This signal contains primary and seconaary 

15 signal portions • The second input to the joint process 

estimator i s th e ~117 -xefexence s ignal vectors -of 27C--3a-mpi-e'3 — 
each. 

The joint process estimator also receives a lambda input 
543, a minimum error input 544 and a number of cells 

20 configuration input 545. These parameters are . .well 

understood in the art. The lambda parameter is often called 
the "forgetting parameter" for a joint process estimator. 
The lambda input 543 provides control for the rate of 
cancellation for the joint process estimator. In the present 

25 embodiment, lambda is set to a low value such as 0.8. 

Because statistics of the signal are non-stationary, a Ixyvr 
value improves tracking. The minimum error input 544 
provides an initialization parameter (conventionally known as 
the "initialization value") for the joint process estimatoo:^ 

30 550* In the present embodiment, the minimum error value isr 

10*^. This initialization parameter prevents the joint 
process estimator 500 from dividing by zero upon initiari 
calculations. The number of cells input 545 to the joint 
process estimator 550 configures the number of cells for the 

35 joint process estimator. In the present embodiment, the 

number of cells for the saturation transform operation 406 is 
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six. As well understood in the art, for each sine wave, the 
joint process estimator requires two cells. If there are two 
sine waves in the 35 - 250 beats/minute range, six cells 
allows for the two heart beat sine waves and one noise sine 
5 wave . 

The joint process estimator 550 subjects the first input 
vector on the first input 542 to a correlation cancellation 
based upon each of the plurality of reference signal vectors 
provided in the second input 540 to the correlation canceler 

10 531 (all 117 reference vectors in sequence in the present- 

embodiment) • The correlation cancellation results in a 
single output vector for each of the 117 reference vectors. 
Each output vector represents the information that the first 
input vector and th-^ ccrr-espcnding re-f-e^rence signal v€rctcr-dc- 

15 not have in common. The resulting output vectors are 

provided as an output to the joint process estimat.ojc. axid- 
subjected to the low pass filter module 552, In the present 
embodiment, the low pass filter 552 comprises a FIR filter 
with 25 taps and with a corner frequency of 10 Hz with the 

20 sampling frequency of 62,5 Hz (i.e., at the decimation 

frequency) . 

The joint process estimator 550 of the present 
embodiment has a settling time of 150 data points. 
Therefore, the last 120 data points from each 270 point 

25 output vector are used for further processing. In the 

present embodiment, the output vectors are further processed 
together as a whole, and are divided into a plurality of bins 
of equal number of data points. As depicted in Figure 18, 
the output vectors are provided to a master power curve 

30 module 554 and to a Divide into five Equal Bins module 556.- 

The Divide into Five Equal Bins module 556 divides each of 
the output vectors into five bins of equal number of data 
points (e.g., with 120 data points per vector, each bin has 
24 data points) . Each bin is then provided to the Bin Power 

35 Curves module 558. 
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The Master Power Curve module 554 performs a saturation 
transform as follows: for each output vector, the sum of the 
squares of the data points is ascertained. This provides a 
sum of squares value corresponding to each output vector 
5 (each output vector corresponding to one of the saturation 

scan values) . These values provide the basis for a master 
power curve 555, as further represented in Figure 22. The 
horizontal axis of the power curve represents the _saJairation 
axis scan values and the vertical axis represents the-^r» of 

10 squares value (or output energy) for each outpttt vector • In- 
other words, as depicted in Figure 22, each of the sum of 
squares could be plotted with the magnitude of the sum of 
squares value plotted on the vertical "energy output" axis aH 
the point on "the* 'hoiizouLal axis of- Llit= -ctrnrreopundxng" 

15 saturation scan value which generated that output vector. 

This results in -a-mas4:er-pQMer-cu3:3/-e ~SS8-,-~a,n~ e^ ?Mi p2j? ^fr- ^hSj::y: ^ 
is depicted in Figure 22. This provides a saturation 
transform in which the spectral content of the attenuated 
energy is examined by looking at every possible saturation 

20 value and examining the output value for the assumed 

saturation value. As will be understood, where the firs-t and 
second inputs to the correlation canceler 531 are mostly 
correlated, the sum of squares for the corresponding output 
vector of the correlation canceler 531 will be very low. 

25 Conversely, where the correlation between the first and 

second inputs to the correlation canceler 531 are not 
significantly correlated, the sum of squares of Lhe^ ouLnnL- 
vector will be high. Accordingly, where the spetrtrrarl ownLcuL 
of the reference signal and the first input to the. 

30 correlation canceler are made up mostly of physiological 

(e.g., movement of venous blood due to respiration) and non- 
physiological (e.g. / motion induced) noise, the output -ener-gy 
will be low. Where the spectral content of the reference 
signal and the first input to the correlation canceler are 

35 not correlated, the output energy will be much higher. 
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A corresponding transform is completed by the Bin Power 
Curves module 558, except a saturation transform power curve 
is generated for each bin- The resulting power curves are 
provided as the outputs of the saturation transform module 
406. 

Tn general, in accordance with the signal model of the 
present invention, there will be two peaks in the power 
curves, as depicted in Figure 22. One peak corresponds to 
the arterial oxygen saturation of the blood, and one peak 
corresponds to the venous oxygen concentration of the bloods. 
With reference to zhe signal model of the present inve nlri o n; 
the peak corresponding to the highest saturation value (noc 
necessarily the peak with the greatest magnitude) correspoirdsr 

to the pr^^ p^^^"^ coefficient r^* Tn nrhp-*" words, the. 

proportionality coefficient r^ corresponds to che- 
red/infrared ratio which will be measured for the arterial 
saturation. Similarly, peak that corresponds to the lowest 
saturation value (not necessarily the peak with the lowest 
magnitude) will generally correspond to the venous oxygen 
saturation, which corresponds to the proportionality 
coefficient r^ in the signal model of the present invention. 
Therefore, the proportionality coei*fdcxeiit — r;; — wi-11 — be — a 
red/infrared ratio corresponding to the venous oxygen 
saturation. 

In order to obtain arterial oxygen saturation, the peak 
in the power curves corresponding to the highest saturation 
value could be selected. However, to improve confidence in 
the value, further processing is comp1 pred .. E±gur.e,-A9 
illustrates the operation of the saturation calculation 
module 408 based upon the output of the saturation transform 
module 406 and the output of the statistics module 404. As 
depicted in figure 19, the bin power curves and the bin 
statistics are provided to the saturation calculation module 
408, In the present embodiment, the master power curves-are 
not provided to the saturation module 408 but can be 
displayed for a visual check* on system operation. The bin 
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statistics contain the red and infrared RMS values, the seed 
saturation value, and a value representing the cross - 
correlation between the red and infrared signals from the 
statistics module 404. 

The saturation calculation module 408 first determines 
a plurality of bin attributes as represented by the Compute 
Bin Attributes module 560. The Compute Bin Attributes module 
560 collects a data bin from the information from the bin 
power curves and the information from the bin statistics, Ir. 
the present embodiment, this operation involves- placing- irhe- 
saturation value of the peak from each power curve 
corresponding to the highest saturation value in the data 
bin. In the present embodiment, the selection of the "highest 
peak is performed -by f i-rst computing the first •d:e^iv=etive-=v>^ 
the power curve in question by convolving the power curve 
with a <=rr.r.r.i--h^-nrr ^ -i f f >»r■^■n^ i a tor f 1 1 fpr fuTirrinn . Tn the 
present embodiment, the smoothing differentiator filter 
function (using a FIR filter) has the following coefficients: 
0.014964670230367 
0.098294046682706 
0.204468276324813 
2.717182664241813 
5,704485606695227 
0 .000000000000000 
-5.704482606695227 
-2.717182664241813 
-0 .204468276324813 
-0.098294046682706 
-0.014964670230357 
This filter performs the differentiation and smoothing. 
Next, each point in the original power curve in question is 
evaluated and determined to be a possible peak if the 
following conditions are met: (1) the point is at least 2% of 
the maximum value in the power curve; (2) the value of titer 
first derivative changes from greater than zero to less than 
or equal to zero. For each point that is found to be a 
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possible peak, the neighboring points are examined and the 
largest of the three points is considered to be the true 
peak. 

The peak width for these selected peaks is also 
calculated. The peak width of a power curve in question is 
computed by summing all the points in the power curve and 
subtracting the product of the minimum value in the power 
curve and the number of points in the power .curve-. In the- 
present embodiment, the peak width calculation is applied to 
each of the bin power curves. The maximum value is selected 
as the peak width. 

In addition, the infrared RMS value from the entire 
snapshot, the red RMS value, the seed saturaHion value ^rm 
each bin, and the cross correlation between the red and 
infrared signals from the statistics module 4r04 are also 
placed in the data bin. The attributes 3xe then used -to 
determine whether the data bin consists of acceptable data, 
as represented in a Bin Qualifying Logic module 562. 

If the correlation between the red and infrared signals 
is too low, the bin is discarded. If the saturation value of 
the selected peak for a given bin is lower than the seed 
saturation for the same bin, the peak is repl a rpr i wirh fhp 
seed saturation value. If either red or infrared RMS value 
is below a very small threshold, the bins are all discarded, 
and no saturation value is' provided, because the measured 
signals are considered to be too small to obtain meaningful 
data. If no bins contain acceptable data, the exception 
handling module 563 provides a message to the display 33 6 
that the data is erroneous. 

If some bins qualify, those bins that qualify as having 
acceptable data are selected, and those that do not qualify 
are replaced with the average of. the bins that are accepted. 
Each bin is given a time stamp in order to maintain the time 
sequence. A voter operation 565 examines each of the bins 
and selects the three highest saturation " values . These 
values are forwarded to a clip and smooth operation 566. 
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The clip and smooth operation 566 basically performs 
averaging with a low pass filter. The low pass filter 
provides adjustable smoothing as selected by a Select 
Smoothing Filter module 568, The Select Smoothing Filter 
5 module 568 performs its operation based upon a confidence 

determination performed by a High Confidence Test module 570, 
The high confidence test is an examination of the peak width 
for the bin power curves. The width of the peaks provides 
some indication of motion by the patient wider pearks- 

10 indicating motion. Therefore, if the peaks are wide, tlie 

smoothing filter is slowed down. If peaks are narrow, the 
smoothing filter speed is increased. Accordingly, the 
smoothing filter 566 is adjusted based on the confidence 
level. The oucpuc of che clip and smooth module 566 provides 

15 the oxygen saturation values in accordance with the present 

invention. 

In the presently preferred embodiment, the clip and 
smooth filter 566 takes each new saturation value and 
compares it to the current saturation value. If the 

20 magnitude of the difference is less than 16 (percent oxygen 

saturation) then the value is pass. Otherwise, if the ne^v 
saturation value is less than the filtered saturation value, 
the new saturation value is changed to 16 less than the 
filtered saturation value. If the new saturation value is 

25 greater than the filtered saturation value, then the new 
saturation value is changed to 16 more than the filtered 
saturation value. 

During high confidence (no motion) , the smoothing f ilLyr 
is a simple one-pole or exponential smoothing filter which is 

30 computed as follows: 

y(n) - 0.6 * x(n) + 0.4 * y(n-l) 
where x(n) is the clipped new saturation value, and y(n) is 
the filtered saturation value. 

During motion condition, a three-pole IIR (infinite 

35 impulse response) filter is used. Its characteristics are 

controlled by three time constants t^, t^, and t^ with values 
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of 0,985, 0.900, and 0.94 respectively. The coefficients for 
a direct form I, IIR filter are computed from these time 
constants using the following relationships: 
ao = 0 

5 = tb + (tc) (t^ + tb> 

a^ = i-t^) (tj (t^ + tb + (tj (tj) 
a3 = (tb)^(tjMtJ 

bo = 1 - tb - (t^) (t, + (tj (tb) ) 
10 bi = 2 (tb) (tj (t^ - 1) 

b^ - (tb) (te) (t^ -K (t^) (tj . (tb) (t^) (t.) - tj 

Figure 20 and 21 illustrate the pulse rate module 41C 
(Figure 14) in greater detail. As illustrated in Figure 20, 

15 the heart rate module 410 has a transient removal and- 

bandpass filter module 578, a motion artifact suppression 
module 580, a saturation equation module 582, a motion status 
module 584, first and second spectral estimation modules 586, 
588, a spectrum analysis module 590, a slew rate limiting 

20 module 592, an output filter 594, and an output filter 

coefficient module 596. 

As further depicted in Figure 20, the heart r-a-t-e--modul-c- 
410 accepts the infrared and red 570-sample snapshots from 
the output of the decimation module 402. The heart rate 

25 module 410 further accepts the saturation value which is 

output from the saturation calculation module 408. In 
addition, the maximum peak width as calculated by the 
confidence test module 570 (same as peak width calculation 
described above) is also provided as an input' to the heart 

30 rate module 410. The infrared and red sample packets, the 

saturation value and the output of the motion status module 
584 are provided to the motion artifact suppression module 
580 . 

The average peak width value provides an input to a 
35 motion status module 584. In the present embodiment, if the 

peaks are wide, this is taken as an indication of motion. If 
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motion is not detected, spectral estimation on the signals is 
carried out directly without motion artifact suppression. 

In the case of motion, motion artifacts are suppressed 
using the motion artifact suppression module 580. The motion 
5 artifact suppression module 580 is nearly identical to the 

saturation transform module 406. The motion artifact 
suppression module 580 provides an output which connects as 
an input to the second spectral estimation module 588. The 
first and second spectral estimation modules SB£-, 588 hav^^ . 

10 outputs which provide inputs to the spectrum analy^sis -module- 
590. The spectrum analysis module 590 also receives an inpuc 
which is the output of the motion status module 584. The 
output of the spectrum analysis module 590 is the initrrai 
heart rate r^f^^ ^'»*Tn^ na^ <nn of the heart xbJzs. m/^/^iO ^ ^ ^ and— 

15 provided as input to the slew rate limiting -module 532. The 

slew rate limiting module 592 connects to the output filter 
594. The output filter 594 also receives an input from the 
output filter coefficient module 595. The output filter 594 
provides the filtered heart rate for the display 33 6 (Figure 

20 11) , 

In the case of no motion, one of the signals (the 
infrared signal in the present embodiment) is subjecced i:o jjC 
removal and bandpass filtering as represented in the DC 
removal and bandpass filter module 578. The DC removal and 

25 bandpass filter module 578 provide the same filtering as the 
DC removal and bandpass filter modules 536, 538. During no 
motion conditions, the filtered infrared signal is provided 
to the first spectral estimation module 586. 

In the present embodiment, the spectral esti nration 

30 comprises a Chirp Z transform that provides a f rpqiienr y 

spectrum of heart rate information. The Chirp Z transform is 
used rather than a conventional Fourier Transform bprausp. a. 
frequency range for the desired output can be designated in 
a Chirp Z transform. Accordingly, in the present embodiment, 

3 5 a frequency spectrum of the heart rate is provided between 3 0 

and 250 beats/minute. In the present embodiment, the 
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frequency spectrum is provided to a spectrum analysis module 
590 which selects the first harmonic from the spectrum as the 
pulse rate. Usually, the first harmonic is the peak in the 
frequency spectrum that has the greatest magnitude and 
5 represents the pulse rate. However, in certain conditions, 

the second or third harmonic can exhibit the greater 
magnitude. With this understanding, in order to select the 
first harmoniCi. the first peak that has an amplitude of at 
least l/20th of the largest peak in the -spectrum is- 
10 selected) . This minimizes the possibility of select iiivj -db- 

the heart rate a peak in the Chirp Z transform caused by 
noise . 

In the case of motion, a motion artifact suppression is 
completed on che snapshot witii the tuotlwiV a rtif ac t — 

15 suppression module 580, The motion artifact suppression 

module 580 i-s-dep-ict-ed in greater det^i^--i2a--Eig\.ire ?"> - B^-... 

can be seen in Figure 21, the motion artifact suppression 
module 580 is nearly identical to the saturation transform 
module 406 (Figure 18) . Accordingly, the motion artifact 

20 suppression module has a motion artifact reference processor 

570 and a motion artifact correlation canceler 57.1-^ 

The motion artifact reference processor 570 is the same 
as the reference processor 53 0 of the saturation transform 
module 406. However, the reference processor 570 utilizes 

25 the saturation value from the saturation module 408, rather 

than completing an entire saturation transform with the 117 
saturation scan values. The reference processor 570, 
therefore, has a saturation equation module SBl,^ a- r ef-crrefrear-- 
generator 582, a DC removal module 583, and a bandpass fiit-^ 

30 module 585. These modules are the same as corresponding,^ 

modules in the saturation transform reference processor 530. 
In the present embodiment, the saturation equation modul-e-5Sl 
receives the arterial saturation value from the saturation 
calculation module 408 rather than doing a saturation ax±s~ 

35 scan as in the saturation transform module 406. This is 

because the arterial saturation has been selected, and there 
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is no need to perform an axis scan. Accordingly, the output 
of the saturation equation module 581 corresponds to the 
proportionality constant r^ (i.e., the expected red to 
infrared ratio for the arterial saturation value) . 
5 Otherwise, the reference processor 570 performs the same 
function as the reference processor 530 of the saturation 
transform module 406. 

The motion artifact correlation canceler 571 is also 
similar to the saturation transform correlation canceler 5TI 

10 (Figure 18) . However, the motion artifact suppression 

correlation canceler 571 uses a slightly different motion 
artifact joint process estimator 572, Accordingly, the 
motion artifact suppression correlation canceler 571 has a 
joint process estimator 572 and a low-pass filter 5T3r. "Tne 

15 motion artifact joint process estimator 572 differs from the 

saturation transform — JuluL — procegs -estdmatx^x — 55rC- -in - that-- 
there are a different number of cells (between 6 and 10 in 
the present embodiment) , as selected by the Number of Cells 
input 574, in that the forgetting parameter differs (.98 in 

20 the present embodiment) , and in that the time delay du"e "to 

adaptation differs. The low-pass filter 573 is the same— a-s 
the low pass filter 552 of the saturation transform 
correlation canceler 531, 

Because only one saturation value is provided to the 

25 reference processor, only one output vector of 270 samples 

results at the output of the motion artifact suppressiorr 
correlation canceler 571 for each input packet of 570 
samples. In the present embodiment, where ' the infrared 
wavelength is provided as a first input to the correlation 

30 canceler, the output of the correlation canceler 571 provides" 

a clean infrared waveform. It should be understood that, as 
described above, the infrared and red wavelength signalsr 
could be switched such that a clean red waveform is^ provided 
at the output of ' the motion artifact suppression correlation 

35 canceler 571. The output of the correlation canceler 571 is 

a clean waveform because the actual saturation value of the 
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patient is known which allows the reference processor 570 to 
generate a noise reference for inputting to the correlation 
canceler 571 as the reference signal. The clean waveform at 
the output of the motion artifact suppression module 580 is 
5 a clean plethysmograph waveform which can be forwarded to the 

display 336. 

As described above, an alternative joint process 
estimator uses the QRD least squares lattice approach (Figure 
8a, 9a and 10a) . Accordingly, the joint process estimatoj^ 

10 573 (as well as the joint process estimator 550) could be- 

replaced with a joint process estimator executing the QRD 
least squares lattice operation - 

Figure 21a depicts an alternative embodiment of the 
motion artifact Gup?r-eas-icn racdule with a joint procssc- 

15 estimator 572a replacing the joint process estimator 572. 

The joint process estimator 572a comprises a .QRD least 
squares lattice system as in Figure 10a. In accordance with 
this embodiment, different initialization parameters are used 
as necessary for the QRD algorithm. 

20 The initialization parameters are referenced in Figure 

21a as "Number of Cells," "Lambda," "MinSumErr, " "Gamslnit," 
and "SumErrlnit . " Number of Cells and Lambda correspond rcr 
like parameters in the joint process estimator 572. Gamslnit 
corresponds to the 7 initialization variable for all stages 

25 except the zero order stage, which as set forth in the QRD 
equations above is initialized to '1'. SummErrlnit provides 
the 6 initialization parameter referenced above in the QRD 
equations. In order to avoid overflow, the larger of the 
actual calculated denominator in each division in the QRD. 

30 equations and MinSumErr is used. In the present embodiment, 

the preferred initialization parameters are as follows: 
Number of Cells = 6 
Lambda = . 8 
MinSumErr = 10'^*' 

35 Gamslnit = 10'^ 

SumErrlnit = 10'*. 
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The. clean waveform output from the motion artifact 
suppression module 580 also provides an input to the second 
spectral estimation module 588, The second spectral 
estimation module 588 performs the same Chirp 2 transform as 
5 the first spectral estimation module 586. In the case of no 

motion, the output from the first spectral estimation module 
586 is provided to the spectrum analysis module 586; in the 
case of motion, the output from the spcnnd .<=;nprrrr?1 
estimation module 588 is provided to a spectrrtrm: ^Sv^alysd-w 

10 module 590, The spectrum analysis module 590 examines the ■ 

frequency spectrum from the appropriate spectral estimation 
module to determine the pulse rate. In the case of motion/ 
the spectrum analysis module 590 selects the peak ±a the 
spectrum with the highesr amplitude, faecHtrse Irhc mvjulwii 

15 artifact suppression module 580 attenuates all other 

frequencies to -a-vciue-bslcw-the -actual heart rate -peak. ^Itl 
the case of no motion, the spectrum analysis module selects 
the first harmonic in the spectrum as the heart rate as 
described above. 

20 The output of the spectrum analysis module 590 provides 

the raw heart rate as an input to the slew rate limiting 
module 592, which provides an input to an output filter 594. 
In the present embodiment, the slew rate limiting module 592 
prevents changes greater that 20 beats/minute per 2 second 

25 interval. 

The output filter 594 comprises an exponential smoothing 
filter similar to the exponential smoothing filter described 
above with respect to the clip and smooth f xLLci . TItct 
output filter is controlled via an output filter coefTici^rrf 

30 module 596, If motion is large, this filter is slowed down, 

if there is little or no motion, this filter can sample much 
faster and still maintain a clean value. The output from the 
output filter 594 is the pulse of the patient, which is 
advantageously provided to the display 33 6. 
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Alternative To Saturation T ransform Module -- Bank 0£ Filters 
An alternative to the saturation transform of the 
saturation transform module 406 can be implemented with a 
bank of filters as depicted in Figure 23. As seen in Figure 
5 23, two banks of filters, a first filter bank 600 and a 

second filter bank 602 are provided. The first filter bank 
600 receives a first measured signal Sxb(t) (the infrared 
signal samples in the present embodiment) on a corresponding 
first filter bank input 604, and the second filter h a pV cno 
10 receives a second measured signal Sx,(t) (the red samples-in- 
the present embodiment) on a corresponding second filter bank 
input 606. In a preferred embodiment, the first and second 
filter banks utilize static recursive polyphase bandpas's- 
f ilters with fixed cerLter frequencies and corner f r.eqiier:cl£C-- 
15 Recursive polyphase filters are described in an article 
Harris, et. al. "Di3.ital Signal Processing, With Efficient 
Polyphase Recursive All-Pass filters" attached hereto as 
Appendix A. However, adaptive implementations are also 
possible. In the present implementation, the recursive 
20 polyphase bandpass filter elements are each designed -to 

include a specific center frequency and bandwidth. 

There are TI filter elements in each filter bank. Each 
of the filter elements in the first filter bank 600 have a 
matching (i.e., same center frequency and bandwidth) filter 
25 element in the second filter bank 602. The center 

frequencies and the corner frequencies of N elements are eaxJa- 
designed to occupy N frequency ranges, 0 to F^, F^ - Fj, F^ - 
Fj, F3 - F4 . . . Fr-1 - Fn as shown in Figure 22,.. 

It should be understood that the number of filirerr 
3 0 elements can range from 1 to infinity. However, in the- 

present embodiment, there are approximately 120 separate 
filter elements with center frequencies spread estenly acrxi^s- 
a frequency range of 25 beats/minute - 250 beats/minute. 

The outputs of the filters contain .information about t-he 
3 5 primary and secondary signals for the first and second 

measured signals (red and infrared in the present example) at 
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the specified frequencies. The outputs for each pair of 
matching filters (one in the first filter bank 600 and one in 
the second filter bank 602) are provided to saturation 
determination modules 610. Figure 23 depicts only one 
saturation determination module 610 for ease of illustration. 
However, a saturation determination module can be provided 
for each matched pair of filter elements for parallel 
processing. Each satnixatrion deteriiiiuaLiuii luodule Lab a ratio- 
module 616 and a saturation equation module 618. 

The ratio module 616 forms a ratio of the second output 
to the first output. For instance, in the present example-, 
a ratio of eHch red- RMS value to each corresponding infrared 
RMS value (Red/IR) is completed in the ratio module 616, 
output of the ratio module 616 provides an input to the 
saturation equation module 618 which references 3. 
corresponding saturation value for the input ratio. 

The output of the saturation equation modules 618 are 
collected (as represented in the histogram module 620) for 
each of the matched filter pairs. However, the data 
collected is initially a function of frequency and 
saturation. In order to form a saturation transform curve 
similar to the cur^-e depicted in Figure 22, -a hi-stogxaa-^ 
the like is generated as in Figure 24. The horizontal axis- 
represents the saturation value, and the vertical axis 
represents a summation of the number of points (outputs from 
the saturation equation modules 618) collected at each 
saturation value. In other words, if the output of the 
saturation equation module 618 for ten different matched 
filter pairs indicates a saturation value di 98%, then a 
point in the histogram of Figure 24 would reflect a value of 
10 at 98% saturation. This results in a curve similar to the 
saturation transform curve of Figure 22. This operation is 
completed in the histogram module 620. 

The results of the histogram provide a power curve 
similar to the power curve of Figure 22. Accordingly, the 
arterial saturation can be calculated from the histogram by 



selecting the peak (greatest number of occurrences in the 
area of interest) - corresponding to the highest saturation 
value (e,g,, the peak /c' in Figure peaks corresponding to 
the highest saturation value peak. Similarly, the venous or 
background saturation can be determined from the histogram by- 
selecting the peak corresponding to the lowest saturation 
value (e.g./ the peak 'd' in Figure 24), in a manner similar 
to the processing in the saturation calculation module 408, 
It should be understood that as an ^^r^r^t-^M^. t-h^ 
histogram, the output saturation (not necessarily a peak in 
the histogram) corresponding to the highest saturation vaTue 
could be selected as the arterial saturation with the 
corresponding ratio representing r^. Similarly, the o uupuc 
saturation correspondirxg to th-e lowes-t sarturaticn value ccv.ld 
be selected as the venous or background saturation with the 
corresponding xat±Q_xapnesenting r... For example, in this 
embodiment, the entry 'a' in the histogram of Figure 24 would 
be chosen as the arterial saturation and the entry in the 
histogram 'b' with the lowest saturation value would be 
chosen as the venous or background saturation. 

ALTERNATIVE DETERMINATION OF COEFFICIENTS r, AND r„ 
As explained above, in accordance with the preaenu 
invention, primary and secondary signal portions, 
particularly for pulse oximetry, can be modeled as follows: 

Sred = Si + Hi (red) (89) 

Sir = Sj 4- n^ (infrared) (90) 

Si - r^Sj and n^ = r^n2 ^ (SX) 

Substituting Equation (91) into Equation (89) provides the 
following : 

Sred = ^^^2 + ^^T^2 (red) (92) 
Note that Sj,ed ^ir used in the model of eqixaJLions 
(89) - (92) . This is because the discussion below is 
particularly directed to blood oximetry. s^^^ and Sj^ 
correspond to S^ and S2 in the preceding text, and the 
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discussion that follows could be generalized for any measure 

signal Sj. and Sj. 

As explained above, determining r. and r^ (which 
correspond to arterial and venous blood oxygen saturation via 
a saturation equation) can be accomplished using the 
saturation transform described above doing a scan of many 
possible coefficients. Another method to obtain r, -and r^ 
based on red and infrared data is to look f or -r. and i.^ which- 
minimize the correlation between s^ and n^, assuming s^ is at 
least somewhat (and preferably sxibstantially) uncorrelated 
with n,, (where k = 1 or 2) . These values can be found by 
minimizing the following statistical calculation function for 
k = 2: 



Correlation I s.,.n.) = lY^ sAS^. S,.^ r^. rJ,xi.jS^.S,.^^.jr^^.^rU 

i (93) 



where i represents time. 

It should be understood that other correlation functions 
such as a normalized correlation could also be used. 

Minimizing this quantity often provides a unique pair -of- 
r, and r^ if the noise component is uncorrelated to the- 
desired signal component. Minimizing this quantity can be 
accomplished by solving Equations (90) and (92) for and nj, 
and finding the minimum of the correlation for possible 
values of r, and r^. Solving for Sj and provides the 
following: 
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inverting the two-by-two matrix provides: 
Thus^ 



1 1 
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or: 
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Preferably, the correlation of equation (93) is enhanced with 
5 a user specified window function as follows: 

C93a) 

The Blackman Window is the presently preferred 
embodiment. It should be understood that there are many 
10 additional functions which minimize the correlation between 

signal and noise* The function above is simply one. Thus, 

Correlation {s^, n^) = |Z I- ^-^^^r^d, ' ^v^ir) S'^redr^^^^^/ » 



If* 



{93b) 
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In order to implement the minimization on a plurality of 
discrete data points, the sum of the squares of the red 
sample points, the sum of the squares of the infrared sample 
points, and the sum of the product of the red times the 
5 infrared sample points are first calculated (including the 
window function, w^) : 

i-l ' 

i-l ' 

IRR = t iS^^ ) w, 



These values are used in the correlation equation {93b) . 

10 Thus, the correlation equation becomes an equation in terms 
of two variables, r, and r^- To obtain r. and r^, an 
exhaustive scan is executed for a good cross-section 'of 
possible values for r^ and r^ (e.g., 20-50 values each 
corresponding to saturation values ranging jfxom 30 - 105) . 

15 The minimum of the correlation function is then selected and 

the values of r, and r^ which resulted in the minimum are 

chosen as r, and r^. 

Once r. and r^ have been obtained, arterial oxygen 
saturation and venous oxygen saturation can be detexmirred-by 
20 provided r^ and r^ to a saturation equation.,' -such as the 

saturation equation 502 of the statistics module 404 wh±ch- 
provides an oxygen saturation value corresponding to tJ^e^ 

ratios r, and r^ 

In a further implementation to obtain r, and r^, the same 
25 signal model set forth above is again used. In order to 

determine r, and r^ in accordance with this implementation, 
the energy in the signal S2 is maximized under the constraint 
that S2 is uncorrelated with n^. Again, this implementation 
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is based upon minimizing the correlation between s and n and 
on the signal model of the present invention where the signal 
s relates to the arterial pulse and the signal n is the noise 
(containing information on the venous blood, as well as 
motion artifacts and other noise) ; r. is the ratio (RED/IR) 
related to arterial saturation and r^ is the ratio (RED/IR) 
related to venous saturation. Accordingly, in this 
implementation of the present invention, r^ and r^ are 
determined such that the energy of the signal s, is maximized 
where s, and nj are uncorrelated. The energy of-the signal 
Sj is given by the following equation: 



(r-rJ2 i-i i-i ' ' i-i 



= [J?,-2r^^+r,=^i?^l (96) 

where Rj is the energy of the red signal, Rj is the energy of 
the infrared signal and Ri,2 is the correlation between the 
red and infrared signals. 

The correlation between and n^ is given by 

1 " \ 
CorrelatiorAs^,!!^) =- r-^ I (S^^ -r^^ ) (— S^ej-^r^Sj.^ ) 

= -r-^—r-, ^-Rx^^r,*rjR,^,-r^,R,] (97) 
(r,-rj2 

As explained above, the constraint is that s^^ and n^ (k 
_ 2 for the present example) are uncorrelated. This 



"decorrelation constraint" is obtained by setting the 
correlation of Equation (97) to zero as follows: 

-R, * (r,+r,) R,, - r,r^ = 0 (98) 

In other words, the goal is to" maximize equation (94) under 
the constraint of equation (98) . 

In order to obtain the goal, a cost function is defined 
(e.g., a Lagrangian optimization in the present -embodiinent ) 
as follows: 

(99) 



where pL is the Lagrange multiplier. Finding the value of r,, 
r^ and that solve the cost function can be accomplished 
using a constrained optimization method such as described .in 
Luenberger, T.inear S- Nonlinear Programming , Addison- Wesley, 
2d Ed. , 1984 . 

Along the same lines, if we assume that the red and 
infrared signals S„a and S„ are non-static, the functions Ri, 
and Ri2 defined above are time dependent. Accordingly, 
with two equations, two unknowns can be obtained by 
expressing the decorrelation constraint set forth in equation 
(98) at two different times. The decorrelartipn corrstraint 
can be expressed at two different times, tj and as 
follows : 

-i?,(t,) + (r,+rjf2,,{t,) -r,r^(t,) =0 (100) 



-i?,(t,) * (r,+rJR,,(t2) - r,r^,(t,) =0 



(101) 



5 



Because equations (100) and (101) are non-linear in r. and r„ 
a change of variables allows the use of linear techniques to 
solve these two equations. Accordingly, with x = r, + r^; 
y = r.rv equations (100) and (101) become 

R,^{t,)x - J2,(t,)y = i?i(tx) (102) 

These equation (102) and (103) can be solved for x .and. y. 
Then, solving for r, and r^ from the changes- tjf ^rariabies- 
equations provides the following: 

r + Jl = X => r3 - xr,+y = 0 (104) 



10 



Solving equation (104) results in two values for r^. In the 
present embodiment, the r^ value that results in x2-r^y>0 is 
selected. If both values of r^ result in x2-r,y>0, the r, 
that maximizes the energy of Sj (Energy (Sj)) at is 
selected, r^ is then substituted into the equations above "to 
15 obtain r.. Alternatively r, can be found directly in the same 

manner r^ was determined. 

Alternative To Saturation Transform - Complex FFT 

The blood oxygen saturation, pulse rate and a clean 
plethysmographic waveform of a patient can also be obtained 

20 using the signal model of the present invention using a 

complex FFT, as explained further with reference to Figures 
25A-25C. In general, by utilizing the signal model of 
equations (89) -(92) with two measured signals, each with a 
first portion and a second portion, where the first portion 

25 represents a desired portion of the signal and the second 

portion represents the undesired portion of the signal, and 
where the measured signals can be correlated with 
coefficients r, and r,, a fast saturation transform on a 
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discrete basis can be used on the sample points from the 
output of the decimation operation 402. 

Figure 25A corresponds generally to Figure 14, with the 
fast saturation transform replacing the previously described 
5 saturation transform. In other words, the operations of 

Figure 25A can replace the operations of Figure 14 • As 
depicted in Figure 25A, the fast saturation transf^orm is 
represented in a fast saturation transform/pulse rate 
calculation module 630. As in Figure 14, the outputs ^-^^ 

10 arterial oxygen saturation, a clean plethysmographic 
waveform, and pulse rate. Figure 25B and 25C illustrate 
additional detail regarding the fast saturation 
transform/pulse rate calculation module 630. As depicted in 
Figure 2SB, the fast saLura tion transform modul-e- -620 hzc- 

15 infrared log and red log modules 64 0, 642 to perform a log 

noTTnal ^7:a^'jr^n-,pq,^r>-^V*g>^^rlfrarpH and r-P^H 1 pg mpHnl f>g 480. 482 
of Figure 17. Similarly, there are infrared DC removal and 
red DC removal modules 644, 646. In addition, there are 
infrared and red high-pass filter modules 64S, 647, window 

20 function modules 648, 640, complex FFT modules 652, 654, 

select modules 653, 655, magnitude modules 656, 658, 
. threshold modules 660, 662, a point -by-point ratio module 
670, a saturation equation module 672, and a select 
saturation module 680. There are also phase modules 690, 

25 692, a phase difference module 694, and a phase threshold 

module 696, The output of the select saturation module G&IL 
provides the arterial saturation on an arterial saturation 
output line 6 82. 

In this alternative embodiment, the snapshot for red and- 

30 infrared signals is 562 samples from the decimation module, 

402. The infrared DC removal module 644 and the red DC 
removal module 646 are slightly different from the infrared- 
and red DC removal modules 4 84, 486 of Figure 17. In the 
infrared and red DC removal modules 644, 646 of Figure 25B, 

35 the mean of all 563 sample points for each respective channel 

is calculated. This mean is then removed from each 



-105- 



individual sample point in the respective snapshot in order 
to remove the baseline DC from each sample. The outputs of 
the infrared and red DC removal modules 644, 646 provide 
inputs to respective infrared high-pass filter module 645 and 
5 red high-pass filter module 647. 

The high-pass filter modules 645, 647 comprise FIR 
filters with 51 taps for coefficients. Preferably, the high- 
pass filters comprise Chebychev filters with a side-lobe 
level parameter of 30 and a corner frequency of .5 Hz (i.e., 
10 30 beats/minute) . It will be understood that this filteji 

could be varied for -performance . With 562 sample points= 
entering the high-panss filters, and with 51 taps for- 
coefficients, there are 512 samples provided from these^ 
respective infrared and red snapshots at the output of the 
15 high-pass filter modules. The output of the high-pass filter 

modules provides an input to the window function modules ff^'S, 
650 for each respective channel. 

The window function modules 648, 650 perform a 
conventional windowing function. A Kaiser windowing function 
20 is used in the present embodiment. The functions throughout 

Figure 25B maintain a point -by-point analysis. In the 
present embodiment, the time bandwidth product for the-K-a-i-se-r 
window function is 7. The output of the window function 
modules provides an input to the respective complex Fast 
25 Fourier Transform (FFT) modules 652, 654. 

The complex FFT modules 652, 654 perform complex FFTs on 
respective infrared and red channels on the data snapshots. 
The data from the complex FFTs is then analyzed in two paths, 
once which examines the magnitude and one which examines the 
30 phase from the complex FFT data points. However, prior to 

further processing, the data is provided to respective 
infrared and red select modules 653, 655 because the output 
of the FFT operation will provide repetitive information from 
0-1/2 the sampling irate and from 1/2 the sampling rate to the 
35 sampling rate. The select modules select only samples from 

0-1/2 the sampling rate (e.g., 0-31.25 Hz in the present 
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embodiment) and then select from those samples to cover a 
frequency range of the heart rate and one or more harmonics 
of the heart rate. In the present embodiment, samples which 
fall in the frequency range of 20 beats per minute to 500 
5 beats per minute are selected. This value can be varied in 
order to obtain harmonics of the heart rate as desired. 
Accordingly, the output of the select modules results in less 
than 256 samples. In the present embodiment:, t:he samp le 
points 2-68 of the outputs of the FFTs are utiXized for 

10 further processing. 

In the first path of processing, the output from the— 
select modules 653, 655 are provided to respective infrared- 
and red magnitude modules 656, 658. The magnitude inodulss^ 
656, 658 perform a magnitude function wherein the macrnitude 

15 on a point-by- point basis of the complex FFT point ? is 

selected for each of the respective channels. The oucpucs of 
the magnitude modules 656, 658 provide an input to infrared 
and red threshold modules 660, 662. 

The threshold modules 660, 662 examine the sample 

20 points, on a point -by-point basis, to select those points 

where the magnitude of an individual point is above a 
particular threshold which is se^t at a pe-rcent-egs- of the 
maximum magnitude detected among all the remaining points in 
the snapshots. In the present embodiment, the percentage for 

25 the threshold operation is selected as 1% of the maximum 

magnitude. 

After thresholding, the data points are forwarded to a 
point -by-point ratio module 670, The point-by -point ratio 
module takes the red over infrared ratio of rhe' values nn a 

30 point-by-point basis. However, a further test is performed 

to qualify the points for which a ratio is taken. As seen in 
Figure 25B, the sample points output from the select modules 
653, 655 are also provided to infrared and red phas^ Tnodules- 
690, 692* The phase modules 690, 692 select the phase value 

3 5 from the complex FFT points. The output of the phase modules 
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690, 692 is then presented to a phase difference module 694. 



The phase difference module 694 calculates the 
difference in phase between the corresponding data points 
from the phase modules 690, 692 • If the magnitude of the. 
phase difference between any two corresponding points is less 
than a particular threshold (e.g., 0.1 radians) in the 
present embodiment) , then the sample points qual±fy. Tf the 
phase of two corresponding sample points is too far apart-, 
then the sample points are not used. The output of the -phase- 
threshold module 695 provides an, enable input to the RED/IR 
rate module 670. Accordingly, in order for the ratio of a 
particular pair of sample points to be taken, the three te s ts 

are executed: 

6 . the red sample must pass the red threshold 

660; 

7. the infrared sample must pass the infrared 

threshold 662; and 

8. the phase between the two points must be less 
than the predefined threshold as determined in the phase 

threshold 696. 

For those sample points which qualify, a ratio is taken 
in the ratio module 670. For those points which do not 
qualify, the saturation is set to zero at the output of the 
saturation equation 672. 

The resulting ratios are provided to a saturation 
equation module which is the same as the saturation equation- 
modules 502, 520 in the statistics module 51^4. In oriTeX" 
words, the saturation equation module 672 accepts the- ratrxrr 
on a point -by-point basis and provides as 3xi output a 
corresponding saturation value corresponding to the discrete 
ratio points. The saturation points outptrtr from the- 
saturation equation module 672 provide a series of saturation 
points which could be plotted aa saturation with respect to 
frequency. The frequency reference was entered into the 
points at the complex FFT stage. 
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The arterial (and the venous) saturation can then be 
selected, as represented in the select arterial saturation 
module 680, in one of two methods according to the present 
invention. According to one method, the arterial saturation 
5 value can be selected simply as the point corresponding to 

the largest saturation value for all points output from the 
saturation equation module 672 for a packet. Alternatively, 
a histogram similar to the histogram of Eigure 22 can be- 
generated in which the number of saturation values -arr" 

10 different frequencies (points) are summed to form a histogram 
of the number of occurrences for each particular saturation 
value • In either method, the arterial saturation can be 
obtained and provided as an output to the select arterial 
saturation module on the arterial saturation output Tine FETZT 

15 In order to nb^a^ n the venous saturation, the minimum 

arterial sacuracion value, o:f" poinLs LhaL eAlilbit iiou -'-geiro" 
value, is selected rather than the maximum arterial 
saturation value- The saturation can be provided to the 
display 336. 

2 0 The fast saturation transform information can also -be 

used to provide the pulse rate and the clean plethysmographic 
wave form as further illustrated in Figure 25C. In order to 
obtain the pulse rate and a clean plethysmographic wave form, 
several additional functions are necessary. As seen in 
25 Figure 25C, the pulse rate and clean plethysmographic wave 

form are determined using a window function module 700, a„ 
spectrum analysis module 702 and an inverse window function 
module 704 , , . . 

As depicted in Figure 25C, the input to the window 

3 0 function module 700 is obtained from the output of che 

complex FFT modules 652 or 654, In the present embodiment^ 
only one measured signal is necessary. Another input to che 
window function module 700 is the arterial saturation 
obtained" from the output of the select arterial saturation 
3 5 module 680. 
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The window function module performs a windowing function 
selected to pass those frequencies that significantly 
correlate to the frequencies which exhibited saturation 
values very close to the arterial saturation value. In the 
present embodiment, the following windowing function is 
selected: 



1- 



100 



where SAT^ equals the saturation value corresponding to each 
particular frequency for the sample points and SAT^ 

10 represents the arterial saturation as chosen at the output of 

the select arterial saturation module 680, This window 
function is applied to the window function input representing 
the complex FFT of either the red or "the "iirrrartid aiyiial-r 
The output of the window function module 700 is a red or 

IS infrared signal represented with a frequency spectrum as 

determined by the FFT, with motion artifacts removed by the 
windowing function. It should be understood that .many 
possible window functions can be provided. In addition, wi-th 
the window function described above, it should he understood 

20 ' that using a higher power will provide more noise 
suppression* 

In order to obtain pulse rate, the output points from 
the window function module 700 are provided to a spectrum, 
analysis module 702. The spectrum analysis module 702 is the^ 

25 same as the spectrum analysis module 590 of Figure 20, In 

other words, the spectrum analysis module /uz determines the 
pulse rate by determining the first harmonic in the frequency 
spectrum represented by the output points of the windowing 
function 700, The output of spectrum analysis module TUT is 

3 0 the pulse rate. 

In order to obtain a clean plethysmographic waveform^ 
the output of the windowing function 700 is applied to an 
inverse window function module 704. The inverse window 
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function module 704 completes an inverse of the Kaiser window 
function of the window function module 648 or 650 of Figure 
25B. In other words, the inverse window function 704 does a 
point -by-point inverse of the Kaiser function for points that 
5 are still defined. The output is a clean plethysmographic 

waveform. 

Accordingly, by using a complex FFT and windowing 
functions, the noise can be suppressed from " the 
plethysmographic waveform in order to obtain the arterial- 

10 saturation, the pulse rate, and a clean pi e t hysmog r aphi c 
waveform. It should be understood that although the above 
description relates to operations primarily in the frequency 
domain, operations that obtain similar results could also be 
accomplished in the time domain. 

15 Relation to Generalized Equations 

The measn re inen^g HAgorih AH for ntn hp oximetry above are 
now related back to the more generalized discussion above. 
The signals (logarithm converted) transmitted through the 
finger 310 at each wavelength Xa and Xb are: 

20 

Sxa ( t ) =Sxredl ( t ) = £10=02. XaC*Hb02X* ( t ) + 6^,. XaC^Kb^'' ( t ) 

+ €„«,,,,,c^«>o2X^(t) +6„,,,.cVx^(t) +n,.{t) ; dCS) 
25 S,,{t)=€„b02.x.c*Hb02X*(t)+£ffl,.xaC*„i,xMt) + n,,(t) ; Om 

Sxa(t)= s,,(t) + n,a(t) ; OTEdj 



30 



35 



Sxb(t:)=£„bo2.xbC\bo2xMt)+e„b.xbC^HbX*(t) +nxb{t) 

Sxb(t)= Sxb(t) + nxb(t:) <^ 
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The variables above are best understood as correlated to 
Figure 6c as follows: assume the layer in figure 6c 
containing Aj and A« represents venous blood in the test 
medium, with A3 representing deoxygenated hemoglobin (Hb) and 
At representing oxygenated hemoglobin (HB02) in the venous 
blood. Similarly, assume that the layer in Figure 6c 
containing A5 and A< represents arterial blood in the test 
medium, with Ag representing deoxygenated hemoglobin (Hb) and 
A< representing oxygenated hemoglobin (HB02) in the arterial 
blood. Accordingly, c''Hb02 represents the concentration of 
oxygenated hemoglobin in the venous blood, c*Hb represents- 
the concentration of deoxygenated hemoglobin in the venous 
blood, x^ represents the thickness of the venous blood (e.g., 
the thickness the layer containing A3 and A,-)- jsimilajiLy r 
c*Hb02 represents the concentration of oxygenated hemoglobin 
in the arterial blood, c*Hb represents the concentration of 
deoxygenated hemoglobin in the arterial blood, and x* 
represents the thickness of the arterial blood (e.g., the 
thickness of the layer containing A5 and A^) 

The wavelengths chosen are typically one in the visible 
red range, i.e., Xa, and one in the infrared range, i.e., Xb. 
Typical wavelength values chosen are Xa = eoOrntr -and- Xb-^ 
910nm. In accordance with the constant saturation method, it 
is assumed that c^oz (t) /cV(t) = constant^ and c''a,02 < t ) /c''„b ( t ) 
= constantj. The oxygen saturation of arterial and venous 
blood changes slowly, if at all, with respect to the sample 
rate, making this a valid assumption. The proportionality 
coefficients for equations (105) and (106) can thBH ±lb 
written as: 



r,(t) = 




(107) 



,c^pc(t) 



Sxa(t) = r,(t)S;,b(t) 



(108a) 
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nx*(t) = r^{t)nxb(t) (108b) 
Sx*(t) ^ r^(t)Sxb(t) (109b) 



In pulse oximetry, it is typically the case that both 
equations (108) and (109) can be satisfied sijnuliianeously. 
Multiplying equation (106) by r^(t) and then -&ubtra:C ting 
10 equation (106) from equation (105), a non-z ero -sexrondary 

reference signal n' (t) is deteinnined by: 

n' (t) = Sx^(t) - r,{t)Sxb(t) Olffe) 

= ^ Hb02. XaC'^HbOl^'' ( ^ ) + C^a^C^'^X^ ( t ) -^Tl^^ ( t ) 

IS 

-r.(t) C€„i^2,xbC^Hb02X^{t)+e«,,;^cVx^(t)+n,^(t)] • (UOa) 

Multiplying equation (106) by rv(t) and then subtracting 
equation (106) from equation (105) , a non-zero primary 
20 reference signal s' (t) is determined by: 

s'(t) = Sx^(t) - r^(t)Sxb(t) (110b) 



= Sxa(t) - r,(t)s.K(t) (111b) 



25 



The constant saturation assumption does noti cause the 
venous contribution to the absorption to be carxeelcd" arl^iig^ 
with the primary signal portions Sxa{t) and S;^i,{t) . Th ii S ; 
frequencies associated with both the low frequency modulated 

30 absorption due to venous absorption when the patient is still 

and the modulated absorption due to venous absorption -;vhen 
the .patient is moving are represented in the secondary 
reference signal n' (t) . Thus, the correlation canceler~or- 
other methods described above remove or derive both 

35 erratically modulated absorption due to venous blood in the 
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finger under motion and the constant low frequency cyclic 
absorption of venous blood. 

To illustrate the operation of the oximeter of Figure 11 
to obtain clean waveform. Figures 26 and 27 depict signals 
5 measured for input to a reference processor of the present 

invention which employs the constant saturation method, i.e, , 
the signals S^^(t) = S^red^t) S^^i^) = S^iaCt) . A first 

segment 26a and 27a of each of the siynals — is relaLlvely 
undisturbed by motion artifact, i.e., the patient did not 

10 move substantially during the time period in which these 
segments were m e asur ed. These segments 26a axKi 27a a re" ^'nus 
generally re presBxiLat ive of the primary plethysmographl-c- 
waveform at each of the measured wavelengths. A se <-x^d 
segment 26b and 27b of each of the signals is affected by 

15 motion artifact, i.e., the patient did move during the time^ 

period in which these segments were measured. "Each of these 
segments 26b and 27b shows large motion induced excursions in 
the measured signal. A third segment 26c and 27c of each of 
the signals is again relatively unaffected by motion artifact 

20 and is thus generally representative of the primary 

plethysmographic waveform at each of the measured 
wavelengths . 

Figure 28 shows the secondary reference signal n' (t) = 
^^xa(^) ~ ^T^^b^t)' determined by a reference processor of 

25 the present invention. Again, the secondary reference signal 

n' (t) is correlated to the secondary signal portions n^a and 
n^b. Thus, a first segment 28a of the secondary reference 
signal n' (t) is generally flat, corresponding, to the fact 
that there is very little motion induced noise In the f i rst 

30 segments 26a and 27a of each signal, A second segment 28b of 

the secondary reference signal n' (t) exhibits large 
excursions, corresponding to the large motion induced 
excursions in each of the measured signals. A third segment 
28c of the noise reference signal n' (t) is generally flat, 

3 5 again corresponding to the lack of motion artifact in the 

third segments 26c and 27c of each measured signal. 
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It should also be understood that a reference processor 
could be utilized in order to obtain the primary reference 
signal s' (t) = s^^-r^s^^^^it) . The primary reference signal 
s' (t) would be generally indicative of the plethysmograph 
5 waveform. 

Figures 29 and 30 show the approximations s'\^(t) and 
s"xb(t) to the primary signals s^^it) and S;^(t) as estimated 
by a correlation canceler using a secondary rftf erence signal 
n' (t) . Note that the scale of Figures 26 through 3-G is- net 

10 the same for each figure to better illustrate changes in each- 
signal. Figures 29 and 30 illustrate the effect of 
correlation cancellation using the secondary reference signal 
n' (t) as determined by the reference processor. Segments 29h 
and 3 0b are noc dominaced by motion iTTdnced~~iiw±s^ "33 — wtwc;- 

15 segments 26b and 27b of the measured signals. Additionally, 

segments 2-9€iv 30a, 29c, -and 3 0c have not -been— suhF?^an^^ al 1 y^., 
changed from the measured signal segments 26a, 27a, 26c, and 
27c where there was no motion induced noise. 

It should be understood that approximation n'\^(t) and 

20 n"xt>(t) to the secondary signals n^^it) and nxb(t) as estimated 

by a correlation canceler using a primary reference si gnal 
s' (t) can also be determined in accordance with the present 
invention. 

METHOD FOR ESTIMATING PRIMARY AND SECONDARY 
25 SIGNAL PORTIONS OF 

MEASURED SIGNALS IN A PULSE OXIMETER 
Implementing the various embodiments of the correl at ion 
canceler described above in software xs- T^ldLi'rei-y 
straightforward given the equations set forth' above, and ther 
30 detailed description above. However, a copy of a comput-er 

program subroutine, written in the C programming language, 
which calculates a primary reference s' (t) using the constant 
saturation method and, using a joint process estimator 572 
which implements a joint process estimator using the 
35 equations (54) - (64) is set forth in Appendix B. This joint 

process estimator estimates a good approximation to the 
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primary signal portions of two measured signals, each having 
a primary portion which is correlated to the primary 
reference signal s' (t) and a secondary portion which is 
correlated to the secondary reference signal n' (t) . This 
5 subroutine is another way to implement the steps illustrated 

in the flowchart of Figure 9 for a monitor particularly 
adapted for pulse oximetry. The two signals are measured at 
two different wavelengths Xa and Xb, where Xa is typically in 
the visible region and Xb is typically -iiL the. in fearer! 
10 region. For example, in one embodiment of the present 

invention, tailored specifically to perform pulse oximetry 
using the constant saturation method, Xa = 660nm and Xb = 
940nm. 

The correspcnde-nce of the prcgras: variables- tc the 
15 variables defined in equations (54) - (64) in the discussion 

of the joint process estimator is as follows: 

A„{t) = nc [m] .Delta 

r£,„(t) = nc[m3 -fref 
20 rb,„(t) = nctm] .bref 

f„(t) = nctm] .ferr 

b„(t) = nc [m] .berr 

0„(t) = nc[m].Fswsqr 

6„{t) = nc[m].Bswsqr 
25 Tmtt.) = nc[m] .Gamma 
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A first portion of the program performs the 
initialization of the registers 90, 92, 96, and 98 and 
intermediate variable values as in the "INITIALIZED 
CORRELATION CANCELER" action block 120. A second portion of 
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the program performs the time updates of the delay element 
variables 110 with the value at the input of each delay 
element variable 110 is stored in the delay element variable 
110 as in the "TIME UPDATE OF LEFT [Z'M ELEMENTS" action 
5 block 130. The calculation of saturation is performed in a 

separate module. Various methods for calculation of the 
oxygen saturation are known to those skilled in the art. One 
such calculation is described in the articles by G.A. Mook, 
et al, and Michael R. Neuman cited above. Once the- 

10 concentration of oxygenated hemoglobin and deoxygenated 
hemoglobin are determined, the value of the saturation is 
determined similarly to equations (72) through (79) wherein 
measurements at times t^ and tj are made at different, yet 
proximate times ovbt' whdch the saturation is relatively 

15 constant. For pulse oximetry, the average saturation at time 

t - (ti+tj) /2 is -then-determined by.;- 

Sat,^,.rui(t) =-_f^fii£l- (112a) 



^Hb.Xj ^Hb.Xb^^^Xj/^^Xb^ (112b) 



Sat it) = (113a) 

^venous ^ ^ y 



^Hb.xr^m>,xb^^^xJ^^xP {113b) 



A third portion of the subroutine calculates the primary 
20 reference or secondary reference, as in the "CALCULATE 

PRIMARY OR SECONDARY REFERENCE (s' (t) or n' (t) ) FOR TWO 
MEASURED SIGNAL SAMPLES" action block 140 for the signals 
Sxa(t) and S;,b<t) using the proportionality constants r^(t) and 
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r^Ct) determined by the constant saturation method as in 
equation (3) , The saturation is calculated in a separate 
subroutine and a value of r^{t) or ry(t) is imported to the 
present subroutine for estimating either the primary portions 
5 Sxa(t) and s^bCt) or the secondary portions n^aCt) and n;^b(t) of 

the composite measured signals Sj,^{t) and S^bCt) . 

A fourth portion of the program performs Z- stage update 
as in the "ZERO STAGE UPDATE" action block 150 where the Z"- 
stage forward prediction error F^it) and Z-stage backward 

10 prediction error h^it) are set equal to the value of th^ 

reference signal n' (t) or s' (t) just calculated. 
Additionally zero-stage values of intermediate variables ti^ 
and i3o(t) (nc [m] .Fswsqr and nc [m] .Bswsqr in the program) are 
calculated for use ln-.sett ing r-f^g-i g^A^-g a^i. go oi£,, axLd_98. 

15 values in the least-squares lattice predictor TO in the 

regression filters 80a and 80b. 

A fifth portion of the program is an iterative loop 
wherein the loop counter, is reset to zero with a maximum 
of m=NC_CELLS, as in the "m=0" action block 160 in Figure 9. 

20 NC_ CELLS is a predetermined maximum value of iterations for 

the loop. A typical value for NC_CELLS is between 6 and 10, 
for example. The conditions of the loop are ^setr buch Lli<x ' c 
the loop iterates a minimum of five times and continues to 
iterate until a test for conversion is met or m-NC_ CELLS. 

25 The test for conversion is whether or not the sum of the 

weighted sum of four prediction errors plus the weighted sum 
of backward prediction errors is less than a small number, 
typically 0.00001 (i.e., 0^{t)+)Sm(t) s 0.00001). 

A sixth portion of the program calculates the forward 

30 and backward reflection coefficient r^.f{t) and rj„,b(^^ 

register 90 and 92 values (nc [m] .fref and nc [m] .bref in the 
program) as in the "ORDER UPDATE m^^-STAGE OF LSL- PREDICTOR" 
action block 170. Then forward and backward prediction 
errors f^(t) and b„(t) (nc[m].ferr and nc[m].berr in the 

35 program) are calculated. Additionally, intermediate 

variables 0^{t), j3m{t) , and7(t) (nc [m] . Fswsqr , nc [m] . Bswsqr, 
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nc [m] . gamma in the program) are calculated. The first cycle 
of the loop uses the value for nc[0] .Fswsqr and nc[0] ,Bswsqr 
calculated in the ZERO STAGE UPDATE portion of the program, 
A seventh portion of the program, still within the loop 
5 begun in the fifth portion of the program, calculates the 

regression coefficient register 96 and 98 values 
^«.x*(t) and K^^xh(^) (nc [m],K_a and nx:-[ni3 JLJolJji the, program) 
in both regression filters, as in the "ORDER IJPDJrrEF m™ S^^£" 
OF REGRESSION FILTER (S) " action block 180. Intermediate 

10 error signals and vari ables e^^xt^it) ^ ^m.xt(^^ ^ P=L^(t)^ ^r^^. 

Pm,xb(t) {ncCm3.err_a and nc[m].err_b, nclm] .rph_a, apd 
nctml.roh_b in the subroutine) are also calculated^ 

The loop iterates until the test for convergence is 
passed- The test for convergence of the joint process 

15 estimator is performed each time the ' loop iterates 

analogously to the "DONE" action block 190 • If the sum ot 
the weighted sums of the forward and backward prediction 
errors 3F„(t) + iSm(t) is less than or equal to 0,00001, the 
loop terminates. Otherwise, sixth and seventh portions of 

20 the program repeat, - - 

The output of the present subroutine is a good 
approximation to the primary signals s"xa(t) and s^^b^^) 
the secondary signals n"xa(t) and n"xb{t) for the set of 
samples S^aCt) and Sxb(t) input to the program. After 

25 approximations to the primary signal portions or the 

secondary signals portions of many sets of measured signal 
samples are estimated by the joint process estimator, a 
compilation of the outputs provides waves which are good 
approximations to the plethysmographic wave or motion 

30 artifact at each wavelength, Xa and Xb. 

It should be understood that the subroutine of Appendix 
B is merely one embodiment which implements the equations 
(54) - (64) . Although implementation of the normalized and 
QRD-LSL equations is also straightforward, a subroutine for 

35 the normalized equations is attached as Appendix C, and a 
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subroutine for the QRD-LSL algorithm is attached as Appendix 
D. 

While one embodiment of a physiological monitor 
incorporating a processor of the present invention for 
5 determining a reference signal for use in a correlation 

canceler, such as an adaptive noise canceler, to remove or 
derive primary and secondary components from a physiological 
measurement has been described in the form of a _pulse 
oximeter, it will be obvious to one skilled in the art~that_ 

10 other types of physiological monitors may also Bmploy — the==- 
above described techniques. 

Furthermore, the signal processing techniques described 
in the present invention may be used to compute the arterial 
and venous blood oxygen saturations ot a phyolwlw^^jloa ^l ^ys^-.-= 

15 on a continuous or nearly continuous time basis. These 

calculationrs -ma-y-be-pe-rfcrmed, regardl^ess— o£— ss?heth*=*r oT^^^ot., 
the physiological system undergoes voluntary motion. 

Furthermore, it will be understood that transformations 
of measured signals other than logarithmic conversion and 

20 determination of a proportionality factor which allows 

removal or derivation of the primary or secondary g nal 
portions for determination of a reference signal are 
possible. Additionally, although the proportionality factor 
r has been described herein as a ratio of a portion of -a 

25 first signal to a portion of a second signal, a similar 

proportionality constant determined as a ratio of a portion 
of a second signal to a portion of a first signal could 
equally well be utilized in the proce ssus -of*^ tire pire&ent- 
invention. In the latter case, a secondary reference signal- 

30 would generally resemble n' (t) = nxb(t) - rnxa^{t) • 

Furthermore, it will be understood that correlation 
cancellation techniques other than joint proces-s esti-ma-t-icn- 
may be used together with the reference signals of the 
present invention. These may include but are not limited-~tx>^ 

3 5 least mean square algorithms, wavelet transforms, spectral 
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estimation techniques, neural networks, Weiner and Kalman 
filters among others. 

One skilled in the art* will realize that many different 
types of physiological monitors may employ the teachings of 
5 the present invention. Other types of physiological monitors 

include, but are in not limited to, electro cardiographs, 
blood pressure monitors, blood constituent monitors (other 
than oxygen saturation) monitors, capnographs, heart rate 
monitors, respiration monitors, or depth. jcvf_ f^7^^<zi-y^^cz^ ^ 

10 monitors. Additionally/ monitors which measure the pressure- 

and quantity of a substance within the "body such as a 
breathalizer, a drug monitor, a cholesterol monitor, a 
glucose monitor, a carbon dioxide monitor, a glucose monicor, 
or a carbon monoxide monitor may also -e?r;plcy the- ^bcvc— 

15 described techniques- 

Furthermore, one skilled in the art will realize that 
the above described techniques of primary or secondary signal 
removal or derivation from a composite signal including both 
primary and secondary components can also be performed on 

20 electrocardiography (EGG) signals which are derived . from 

positions on the body which are close and highly correlated 
to each other. Tt should be understood that a tripoiar 
Laplacian electrode sensor such as that depicted in Figure 31 
which is a modification of a bipolar Laplacian electrode 

25 sensor discussed in the article "Body Surface Laplacian EGG 

Mapping" by Bin He and Richard J. Cohen contained in the- 
journal IEEE Transactions on Biomedical Engineering, Vol. 39, 
No, 11, November 1992 could be used as an EGG censor. It rr.ust 
also be understood that there are a myriad of possible EGG 

30 sensor geometry's that may be used to satisfy the 

requirements of the present invention. The same type of 
sensor could also be used for EEG and EMG measurements. 

Furthermore, one skilled in the art will realize that 
the above described techniques can also be performed on 

3 5 signals made up of reflected energy, rather than transmitted 

energy. One skilled in the art will also realize that a 
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primary or secondairy portion of a measured signal of any type 
of energy, including but not limited to sound energy. X-ray 
energy, gamma ray energy, or light energy can be estimated by 
the techniques described above. Thus, one skilled in the art 
5 will realize' that the techniques of the present invention can 
be applied in such monitors as those using ultrasound where 
a signal is transmitted through a portion of the body and 
reflected back from within the body back through this portion 
of the body* Additionally, monitors such as echo 

10 cardiographs may also utilize the techniqxjtes=- of the present 

invention since they too rely on transmission and reflection. 

While the present invention has been descril^ed in terms 
of a physiological monitor, one skilled in the art will 
realize that the signal processing t:echnicp,^es of 

15 invention can be applied in many areas, including but not 

limited to the j)rocessin_g_ of a physiological signal. The 
present invention may be applied in any situation where a 
signal processor comprising a detector receives a first 
signal which includes a first primary signal portion and a 

20 first secondary signal portion and a second signal which 

includes a second primary signal portion and a second 
secondary signal portion. Thus, the signal processor of triie 
present invention is readily applicable to numerous signal 
processing areas. 
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